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ABSTRACT 


Nonlinear  motion  equations  were  used  to  predict 
vertical  planar  trajectories  of  the  Basic  Pinner  Missile  ob¬ 
served  previously  in  underwater  model  experiments  at 
California  Institute  of  Technology.  Values  of  ^he  hydrodynamic 
coefficients  appearing  in  the  motion  equations  were  obtained 
from  Basic  Pinner  experiments  made  previously  at  CIT,  Bureau 
of  Standards,  and  Davidson  Laboratory.  Coefficients  also  were 
estimated  on  theoretical  grounds.  The  observed  trajectories 
could  not  be  predicted  using  experimental  coefficients  alone. 
However,  theoretical  coefficients  yielded  predictions  that 
agreed  well  with  the  observed  motions.  Computed  trajectories 
were  particularly  sensitive  to  independent  changes  in  tho 
values  of  the  hydrodynamic  static  and  damping  derivative  coef¬ 
ficients,  However,  when  these  values  were  kept  consistent 
with  relations  expressing  their  physical  interdependence,  they 
could  be  varied  over  fairly  wide  ranges  without  appreciably 
affecting  the  predictions. 
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INTRODUCTION 


The  main  purpose  of  this  study  is  to  determine 
whether  observed  underwater  trajectories  of  the  Basic  Pinner 
Missile  can  be  predicted  using  equations  of  motion  and  hydro- 
dynamic  coefficients  obtained  from  experiment  or  theory.  The 
current  *ork  is  part  of  a  larger  program,  sponsored  by  the 
Bureau  of  Naval  Weapons,  to  verify  present  day  trajectory  pre¬ 
diction  methods. 

The  Basic  Pinner  Missile  (Pig.  l)  consists  of  a 
cylindrical  shaft,  pointed  at  the  forward  end,  and  a  set  of 
large,  wedge-shaped,  cruciform  tail  fins.  Although  used  for 
research  in  both  aerodynamics  and  hydrodynamics,  it  is  essen¬ 
tially  prototypal  of  a  supersonic  aerodynamic  missile.  Partly 
because  of  its  blunt  base,  but  primarily  because  of  the  great 
size  of  its  fins,  the  hydrodynamic  reactions  it  develops  in 
underwater  flight  are  extremely  large  compared  with  those  on 
a  typical  hydrodynamic  body  (e.g.,  missiles,  torpedoes, 
submarines)  of  similar  length  and  cross-sectional  area,  flying 
at  the  same  speed  and  attitude.  Consequently,  the  hydro- 
dynamic  derivative  coefficients  for  the  Basic  Pinner  not  only 
are  very  large  compared  with  those  of  more  representative 
underwater  craft,1  but  also  are  influenced  abnormally  by  the 
lift  and  drag  characteristics  of  Its  fin  configuration. 

Measurements  of  the  hydrodynamic  forces  and  moments 
acting  on  Basic  Finner  models  under  various  conditions  have 
been  made  by  Kermcen,2  Kiceniuk,3'4  Heald  and  Adams,5  and 
Savitsky  and  Prowse.®  Extensive  model  trajectory  measurements 
have  been  reported  by  Price.7  Stubstad,8  Parrish,0  and 
Stubstad  and  Waugh10  also  have  published  results  of  trajectory 
tests.  In  addition,  some  heretofore  unpublished  Basic  Finner 
trajectory  data  have  been  made  available  to  Davidson  Laboratory 
by  the  Bureau  of  Naval  Weapons. 
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Some  wcrk  has  been  done  also  In  comparing  measured 
Basic  Firmer  trajectory  data  with  results  obtained  using 
motion  equations.  Price7  and  Stubstad  and  Waugh10  have  em¬ 
ployed  nonlinear  motion  equations  as  models  for  curve-fits  of 
observed  trajectory  data,  to  define  hydrodynamic  force  func¬ 
tions  numerically.  Stubstad0  has  used  the  solution  of  a 
linear  system  of  motion  equations  as  a  basis  for  curve-fitting 
observed  trajectory  data  to  obtain  the  numerical  valueB  of 
three  physical  constants.  Caster11  also  has  curve-fitted 
observed  trajectory  data,  by  specifying  all  but  one  of  the 
hydrodynamic  terms  appearing  in  a  system  of  nonlinear  motion 
equations  and  varying  the  unspecified  term.  Parrish9  has 
compared  observed  trajectory  data  with  results  computed  using 
a  linear  system  of  motion  equations.  Of  all  the  previous 
wcrk,  however,  only  Parrish's  is  an  actual  comparison  of  com¬ 
puted  and  observed  trajectories,  and  this,  being  a  single 
comparison,  is  not  a  sufficient  basis  for  any  definite  con¬ 
clusions  regarding  either  tne  correct  forms  of  the  motion 
equations  or  the  adequacy  of  the  hydrodynamic  data. 

The  model  trajectory  experiments  forming  the  basis 
for  comparison  with  the  trajectories  predicted  herein  include 
those  reported  by  Price, 7  and  some  supplied  directly  to 
Davidson  Laboratory  by  the  Bureau  of  Naval  Weapons.  Th,-  ex¬ 
periments  were  performed  at  the  Controlled  Atmosphere  Launch¬ 
ing  Tank  of  the  California  Institute  of  Technology  (CIT).7 

The  Price  results  are  vertical  planar  trajectories 
in  which  heavier- than-water  models  are  dropped  from  rest. 

Both  straight  (seven  runs)  and  curved  (eight  runs)  vertical 
trajectories  are  included.  The  results  supplied  by  the 
Bureau  of  Naval  Weapons  are  vertical  planar  trajectories  in 
which  heavier- than-water  (eight  runs),  lighter- than-water 
(two  runs),  and  neutrally  buoyant  (one  run)  models  are  launched 
with  nonzero  horizontal  velocities. 


These  results  all  belong  to  the  class  of  vertical 
planar  trajectories.  For  purposes  of  the  present  study,  they 
were  categorized  further  according  to  the  number  of  degrees 
of  freedom  of  the  motion  and  the  type  of  initial  conditions. 
Four  categories  were  established:  high  speed  vertical  planar 
trajectories,  low  speed  vertical  planar  trajectories,  straight 
vertical  trajectories,  and  a  straight  horizontal  trajectory. 

The  first  and  fourth  categories  inclide  the  trajec¬ 
tories  supplied  by  the  Bureau  of  Naval  Weapons;  the  second  and 
third  consist  of  the  results  reported  by  Price.  Differentia¬ 
tion  between  the  data  from  the  two  sources  is  based  on  the 
prevailing  flow  conditions.  In  each  trajectory  reported  by 
Price,  the  model  is  dropped  from  rest.  Hence,  the  motions 
run  in  or,  at  least,  through  the  laminar  flow  regime.  In  the 
runs  supplied  by  BuWeps,  the  model  is  launched  at  high  enough 
velocity  so  that  the  flow  is  turbulent  throughout  virtually 
the  entire  trajectory.  Moreover,  the  angles  of  attack  at¬ 
tained  in  the  high  speed  (BuWeps)  trajectories  are  restricted 
to  the  range  |a|<6.5°,  whereas  attack  angles  as  high  as  a=30° 
are  attained  in  the  low  speed  (Price)  runs. 

Since  actual  motions  of  full-scale  missiles  are  in¬ 
variably  in  the  turbulent  flow  regime,  the  results  of  the 
present  study  of  primary  practical  Importance  are  those  per¬ 
taining  to  the  high  speed  (BuWeps)  trajectories.  Although 
emphasis  apparently  was  about  evenly  divided  between  high  and 
low  speed  trajectories  in  the  phase  of  the  Basic  Firmer  pro¬ 
gram  devoted  to  vertical  planar  trajectory  measurements,  the 
relative  unimportance  of  the  low  speed  motions  appears  to 
have  been  considered  in  the  planning  of  the  hydrodynamic  force 
and  moment  measurements.  These  essentially  were  restricted  to 
Reynolds  number  in  the  range  0.4x10®  <;  Rg  ^  8.5x10.°  This  is 
similar  to  the  ranges  attained  in  the  high  speed  trajectories, 
but  ignores  the  range  0  <;  R0  ^  0.4x10, 6  included  in  each  of 


R-898 
-  3- 


the  low  speed  trajectories.  Similarly,  the  range  of  attack 
angles  covered  in  the  force  and  moment  measurements, 

-13°  <  a  <  20°,  encompasses  those  attained  in  the  high  speed 
trajectories  but  not  those  of  the  low  speed  runs. 

The  reason  for  differentiating  between  the  one 
degree  of  freedom  and  three  degrees  of  freedom  trajectories 
is  that  with  zero  angle  of  attack  ( one  degree  of  freedom 
trajectories),  the  complexity  of  the  hydrodynrmi'’  reactions 
le  very  much  reduced  due  to  the  symmetry  of  the  flow.  Since 
Reynolds  number  and  angle  of  attack  effects  are  strongly 
coupled,  this  distinction  is  particularly  significant  for  the 
case  of  the  low  speed  trajectories,  where  the  ranges  of  flight 
angle  of  attack  are  greatest  and  where  the  available  hydro- 
dynamic  force  and  moment  data  are  most  limited. 

This  research  was  performed  under  the  sponsorship 
of  the  Bureau  of  Naval  Weapons  under  BuWeps  Contract 
NOW  60-034l-d  and  DL  Project  LM  2328. 
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EQUATIONS  OP  MOTION 


The  motion  equations  used  in  the  prediction  of  the 
Basic  Pinner  model  trajectories  were  derived  from  general 
equations  of  motion  for  a  submerged  body  formulated  by 
Strumpf.12  The  general  equations  were  simplified  by  the  fol¬ 
lowing  assumptions: 

1.  The  mass  and  mass  distribution  of  the  body  do  not 
vary  with  time. 

2.  The  body  possesses  both  horizontal  and  vertical 
planes  of  symmetry. 

3.  No  thrust  or  control  forces  or  moments  act  on 
the  body. 

4.  The  center  of  gravity  lies  on  the  longitudinal 
centerline  of  the  body. 

5.  The  products  of  inertia  are  zero. 

6.  The  moments  of  inertia  about  the  y  and  z  axes 
are  equal.* 

7.  The  body  is  constrained  to  move  in  a  vertical 

plane,  l.e.,  v=v= r=r=^=0=O.  ( l) 

8.  The  hydrodynamic  derivatives  X'  ,  X'  ,  Z!,  M! 

are  equal  to  zero.  ww  qq  q  w 

9.  Hydrodynamic  terms  of  order  three  and  higher 

(e.g.,  terms  of  the  form  Z»  w3,  Z'  qw2)  are 
negligible.  www  qww 

Assumptions  1  through  6  are  based  upon  the  physical 
properties  of  the  tested  Basic  Pinner  model  configurations. 
Assumption  I  restricts  the  general  equations  to  the  three 
degrees  of  freedom  (vertical  plane)  cases  considered  in  the 
present  study.  Assumption  8  is  customary  in  analyses  of  this 
nature  and  is  consistent  with  the  results  of  previous 
research. 12’ 13 


*  The  coordinate  system  and  definitions  of  all  symbols  used 
In  this  report  are  presented  tn  Appendix  A. 
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Assumption  9  was  dictated  by  practical  considera¬ 
tions.  Provided  the  hydrodynamic  forces  and  momentB  are 
analytic  functions  of  the  body  velocity  and  acceleration  com¬ 
ponents,  they  can  be  approximated  theoretically  to  any  accu¬ 
racy  by  Taylor  series  expansions  of  sufficiently  high  order. 

In  practical  application,  however,  the  inclusion  of  Taylor 
series  terms  of  order  higher  than  the  first  is  discouraged  by 
the  limitations  of  available  experimental  and  theoretical 
methods  for  determining  the  numerical  values  of  nonlinear 
hydrodynamic  derivatives.  The  usual  practice  in  past  trajec¬ 
tory  prediction  work  has  been  to  employ  linear  equations. 8,8,14 
The  equations  used  in  the  present  study,  however,  represent  a 
compromise.  Although  most  nonlinear  hydrodynamic  effects  are 
neglected,  all  inertial  and  gravity  forces  and  moments  are  in¬ 
cluded.  The  drawbacks  of  linear  equations  are  reduced,  where¬ 
as  the  amount  of  required  hydrodynamic  data  is  kept  to  a  mini¬ 
mum.  Like  linearized  equations,  however,  the  present  motion 
equations  should  not  be  expected  to  apply  to  cases  where  the 
hydrodynamic  forces  ard  moments  are  highly  nonlinear  and  the 
motions  deviate  appreciably  from  the  point,  u,  v,  w,  ...=u,  0, 
0,  ...  (a=0),  about  which  the  Taylor  series  are  expanded.12 

On  the  basis  of  assumptions  1  through  9»  the  general 
longitudinal  force  equation,  vertical  force  equation,  and 
pitching  moment  equation, 12  respectively,  reduce  to 

m( d+qw-x^q2)  =  -(W-B)sln6+  £  AxX’(u2+w2)+  £  Ax^Xti'ii'Xwqwq) ,  (2) 

m(w-qu-xQq)  =  (W-B)cose  +  £  AxZ^uw  +  £  Axt( Z^<H-Z^uq)  ,  (3) 

Iyy<HmxG(qu-*)=  -WxQcos0+  £  A^NU  4+  £  Axt2M^uq+  £AxtM^UW*(  4) 

All  of  the  quantities  appearing  in  these  equations  are  defined 
in  the  Nomenclature  (Appendix  A)  and  thoroughly  discussed  in 
ref.  12,  which  i3  readily  available.  Accordingly,  no  detailed 
discussion  will  be  presented  here.  However,  the  various  terms 

*  The  subsidiary  relation  required  for  solution  of  these  equations  re¬ 
duces  to  q  «  9  (see  Appendix  C). 
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will  be  described  in  a  general  way.  The  terms  on  the  left- 
hand  aide  of  each  equation  are  inertial  terms)  the  familiar 
"^(kV)"  and  "^(  I<u) "  terms  of  Newton's  second  law.  The  first 
terms  on  the  right-hand  side  of  the  equations  are  forces  and 
moments  arising  due  to  the  influence  of  gravity.  All  other 
terms  are  forces  and  moments  of  hydrodynamic  origin.  Note 
that  the  artifice  of  introducing  the  so-called  "virtual  masses 
and  "virtual  moments  of  inertia,"  common  in  the  formulation  of 
motion  equations  for  hydrodynamic  bodies,  is  avoided.  No  dis¬ 
tinction  is  drawn  between  the  hydrodynamic  reactions 

often  accounted  for  in  this  way,1®'15  and  any 
other  hydrodynamia  forces  and  moments. 

Equations  2  through  4  are  the  equations  of  vertical 
planar  motion  employed  in  this  study.  They  were  used  in  the 
predictions  "or  both  high  speed  and  low  speed  vertioal  planar 
trajectories.  However,  their  applicability  to  the  two  dif¬ 
ferent  cases  is  far  from  equal.  Under  the  conditions  prevail¬ 
ing  in  the  important  case  of  the  high  speed  trajectories,  all 
assumptions  made  in  their  derivation  are  reasonable.  In  these 
runs,  the  flight  angles  of  attack  are  restricted  to  a  range 
(|a|<  6.5°)  in  which  the  available  data5  indicate  that  the 
nonlinear  angle  of  attack  dependence  (or,  equivalently, 
w-dependence)  of  the  hydrodynamic  forces  and  moments  is  not 
significant.  No  Basic  Pinner  force  and  moment  data  for  non¬ 
zero  values  of  pitch  angular  velocity  q  are  available,  but  re¬ 
sults  of  previoup  Davidson  Laboratory  studies15'17  indicate 
that  q-w  coupling  and  nonlinear  q  terms  may  reasonably  be 
neglected  for  the  ranges  of  angular  velocity  (q1  <  0.2)  and 
angle  of  attack  attained  in  the  high  speed  runs. 

The  situation  is  vastly  different  for  the  case  of 
the  low  speed  trajectories.  Angles  of  attack  as  great  as  30° 
are  attained  in  some  of  these  flights.  Under  these  conditions 
higher  order  Taylor  series  terms  definitely  are  required  to 


represent  the  hydrodynamic  reactions  accurately.  Hence  the 
present  mathematical  model  for  the  hydrodynamic  forces  and 
moments  is  clearly  unrealistic  under  the  conditions  prevail- 
in?  in  the  low  speed  vertical  planar  trajectories. 

The  extreme  attack  angles  in  these  runs  invariably 
are  attained  in  conjunction  with  Reynolds  numbers  lower  than 
the  minimum  value  (0.4x10°)  covered  in  the  Basic  Pinner  force 
and  moment  measurements.8”*  Hence  no  experimental  data  for 
direct  evaluation  of  the  nonlinear  hydrodynamic  effects  is 
available  unless  it  can  be  assumed  that  such  effects  are  not 
significantly  Reynolds  number  dependent.  But  such  an  assump¬ 
tion  would  be  in  contradiction  to  previous  experience.18* 10 
In  any  event,  no  Basic  Pinner  force  and  moment  data  are  avail¬ 
able  for  | at  |  >  20°,  at  any  Reynolds  number.  Beceuse  of  the 
limited  practical  Importance  of  such  low  Reynolds  number 
motions  of  hydrodynamic  missiles,  applicable  experimental  data 
for  ether  configurations  similar  to  the  Basic  Pinner  are  not 
available  either.  Pinaily,  reliable  methods  for  evaluating 
the  required  quantities  theoretically  have  not  been  developed. 

Since  the  low  speed  vertical  planar  trajectories 
are  of  little  practical  interest  anyway,  it  was  felt  that 
efforts  to  obtain  additional  detailed  hydrodynamic  force  and 
moment  data  for  the  low  Reynolds  number-high  attack  angle 
regime  were  not  warranted.  However,  since  the  trajectory 
measurements  had  been  completed  already  and  the  ’’esults  were 
available,7  it  was  felt  that  an  attempt  should  be  made  to  pre¬ 
dict  them  using  the  present  notion  equations. 

For  the  straight  horizontal  trajectory  prediction, 
eqs.  2  through  4  are  modified  by  the  substitutions 

w  =  w  =  q=  q  =  6=xQ=W-B  =  0  (5) 
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Only  the  longitudinal  foroe  equation  remains  nontrivial*  re¬ 
ducing  to 

mfli  "  I  V +  £  V*  10  (6) 

Similarly*  for  the  straight  vertical  trajectory  predictions* 
the  substitutions 

w-w-q-q-0,  8  -  -  £  ,  (7) 

yield,  simply* 

mu  -  £  Vi11*  +  I  AxWJA  +  (W‘B)  (8) 

Note  that  the  restrictions  to  one  degree  of  freedom 
motions  yield  equations  (eqs.  6  and  8*  respectively*  for 
horizontal  and  vertical  flight)  which*  although  nonlinear* 
nevertheless  are  soluble  by  exact  analytic  methods*  which  the 
three  degrees  of  freedom  equations  (eqs.  4  through  6)  are  not 
(see  Appendix  C) .  This  simplification  in  computation  is 
matched  by  the  reduction  in  the  complexity  of  the  actual  phy¬ 
sical  situation  resulting  from  the  absence  of  angles  of  attack 
or  angular  velocities. 
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HYDRODYNAMIC  COEFFICIENTS 


'There  are  two  major  souroee  of  error  in  predicting 
underwater  trajectories  with  motion  equations.1*  The  first 
is  in  the  functional  representations  given  to  the  external 
forces  and  moments,  particularly  those  of  hydrodynamic  origin, 
acting  on  the  missile  during  flight  (discussed  in  the  preced¬ 
ing  section) .  The  aeoond  is  in  the  numerical  values  used  for 
the  various  quantities  appearing  in  the  motion  equations. 

Since  such  quantities  as  the  density  of  water,  the  length  of 
the  missile,  its  cross-sectional  area,  weight,  displacement, 
moment  of  inertia,  and  CQ  location  presumably  are  known  with 
great  accuracy,  the  crux  of  this  problem  lies  in  the  hydro- 
dynamic  derivative  coefficients.  For  the  Basic  Pinner,  values 
for  these  quantities  are  available  from  the  results  of  several 
experiments, 2"fl  as  well  as  from  calculations  based  on  a  number 
of  existing  theories. 

The  Important  features  of  each  of  the  experiments 
whose  data  were  used  to  obtain  values  for  the  hydrodynamic 
force  and  moment  derivative  coefficients  are  summarized  in 
Table  I. 

The  experimental  methods  customarily  employed  for 
evaluating  hydrodynamic  force  and  moment  derivative  coeffi¬ 
cients  can  be  divided  into  two  general  classes.  The  first 
class  consists  of  "static"  methods,  in  which  forces  and 
moments  acting  on  a  model  are  measured  directly  as  functions 
of  the  linear  and  angular  velocity  and  acceleration  components 
( and  any  other  pertinent  variables) .  Then  the  derivatives  are 
evaluated  from  the  data  according  to  their  definitions,  e.g., 
is  the  slope  of  the  curve  of  Z’  vs  w‘,  with  u*=u^, 
u»=v»=v'=w'  =  p’ =/>'=.. .  =  o,  evaluated  at  the  point  wf=0.  In 
the  second  general  class  of  experimental  methods,  the  forces 
and  moments  are  not  measured.  Instead,  forced  oscillatory 
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motions  of  ths  model  are  induced*  and  suoh  quantities  as  dis¬ 
placement  response  amplitude  and  phase  lag  are  recorded.  The 
hydrodynamic  foroes  and  moments  are  expressed  as  linear  func¬ 
tions  of  the  velooity  and  acceleration  components  of  the  model 
relative  to  the  fluid*  and  combined  with  similar  linear  repre¬ 
sentations  for  the  applied  foroes  and  moments  of  mechanioal 
origin*  in  a  system  of  linear  motion  equations  whose  solutions 
are  used  to  express  the  hydrodynamic  derivatives  in  terms  of 
the  measured  motion  parameters.  The  second  column  of  Table  I 
characterises  each  of  the  Basic  Pinner  experiments  as  "static" 
or  "osoillatory depending  upon  the  class  to  which  it  belongs. 

The  derivative  coefficients  obtained  from  the  data 
of  each  experiment  are  listed  in  the  third  column  of  Table  1. 
The  derivation  of  these  quantities  from  the  available  data 
consisted,  generally*  of  three  stepst  l)  evaluation  of  deriva¬ 
tives  from  plotted  data)  2)  conversion  to  a  standard  nondlmen- 
sionaliz&tion  basis)  and  3)  adjustment  to  the  x*y*z*  coordinate 
system  with  origin  at  the  CB.  The  fourth  column  of  Table  I* 
outlines  these  steps  in  detail  for  each  result. 

The  values  of  each  of  the  hydrodynamic  derivative 
coefficients  determined  as  indicated  in  Table  I  were  plotted 
versus  Reynolds  number.  Figure  2  shows  such  a  plot  for  the 
case  of  the  zero  order  derivative  (the  drag  coefficient), 
the  derivative  for  which  the  Reynolds  number  data  are  most  ex¬ 
tensive.  Scatter  between  the  various  experiments  at  a  given 
Reynolds  number  is  as  great  as  the  variation  with  Reynolds 
number  within  a  given  experiment.  This  was  found  to  be  the 
general  case  when  values  for  a  coefficient  were  available 
from  more  than  one  source.  On  this  basis,  it  was  assumed  that 
the  influence  of  Reynolds  number  on  the  Basic  Finner  hydro- 
dynamic  derivative  coefficients  is  negligible  for  Reynolds 
numbers  in  the  range,  0.4x10®  ^  Re  £  8.5x10,®  the  total  range 
of  the  various  experiments.  This  covers  the  ranges  of 
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Reynolds  number  attained  in  all  of  the  high  speed  vertical 
planar  trajectories. 

Experimentally  derived  values  for  the  Basic  Firmer 
hydrodynamic  derivative  coefficients  are  presented  in 
Table  II.  With  the  exception  of  the  results  obtained  from 
the  data  of  Heald  and  Adams, 8  these  values  represent  averages 
estimated  from  all  the  data  reported  in  each  reference.  The 
Heald  and  Adams  data  have  been  given  added  weight  by  consider¬ 
ing  separately  the  results  of  experiments  using  (four)  dif¬ 
ferent  models.  The  emphasis  on  these  data  is  in  accordance 
with  the  relative  thoroughness  and  clarity  of  the  results  pre- 
sented  in  the  various  references.*”® 

In  general,  there  are  substantial  differences  in 
measurements  of  hydrodynamic  force  and  moment  derivative  co¬ 
efficients  for  a  supposedly  identical  configuration  tested  by 
different  laboratories. *°»ai  The  present  case  is  by  no  means 
an  exception  to  this  rule.  The  ranges  of  valuta  determined 
from  the  various  experiments  for  X£,  Z^,  M^,  and  (the  de¬ 
rivatives  for  which  more  than  one  value  are  available)  are, 
respectively,  53#,  20 #,  31#,  and  133#  of  the  average  value  for 
each.  Perhaps  more  significantly,  the  ranges  of  these  values 
determined  from  the  results  of  a  single  laboratory,  the 
National  Bureau  of  Standards  (reported  by  Heald  and  Adams5), 
are,  respectively,  40#,  19#,  8#,  and  44#  of  these  same  averages. 

There  are  three  major  factors  that  generally  con¬ 
tribute  to  the  discrepancies  between  repeated  experimental  de¬ 
terminations  of  the  hydrodynamic  derivative  coefficients.  All 
undoubtedly  are  significant  in  the  case  of  the  Basic  Firmer 
results.  The  first  is  in  the  subjective  methods  used  in 
handling  the  experimental  data.  Each  of  the  values  listed  in 
Table  II  was  established  on  the  basis  of  at  least  one  curve¬ 
fitting  procedure.  Some  of  these  curve  fits  were  performed 
at  Davidson  Laboratory  in  the  course  of  the  present  study 
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(see  Appendix  B) ;  others  were  performed  at  the  faollltlos 
where  the  experiments  were  made.  In  some  of  the  latter 
caeea,a*a>4  It  is  impossible  to  assess  the  importance  of  this 
factor  directly,  because  the  raw  test  data  and/or  descriptions 
of  the  curve-fitting  techniques  used  are  not  available.  It  is 
reasonable  to  assume,  however,  that  the  uncertainty  in  fitting 
these  data  is  no  less  than  that  in  fitting  the  other  Basic 
Firmer  data.  That  this  uncertainty  is  significant  can  be  seen 
dearly  from  the  sample  derivations  of  values  for  hydrodynamic 
coefficients  from  test  data5  presented  in  Appendix  B.  In  fact, 
listing  ranges  of  values  in  Table  II  might  have  been  more  ap¬ 
propriate  than  showing  the  discrete  values.  Thus,  -11.9 
(Table  II,  row  1,  column  l)  would  be  replaced  by  -12.0£Z^£-11.8, 
Similarly,  -2.4  (row  1,  column  2)  would  be  replaced  by 
-2.5<H^<-2.2.  However,  it  was  felt  that  the  uncertainty  in 
the  tabulated  values  would  be  sufficiently  apparent  without 
introducing  this  additional  complication. 

The  second  factor  contributing  to  the  discrepancies 
between  the  repeated  determinations  of  the  hydrodynamic  coef¬ 
ficients  is  actual  physical  differences  in  the  hydrodynamic 
force  and  moment  components  measured  under  supposedly  identi¬ 
cal  conditions.  Such  differences  can  arise  from  many  sources. 
Perhaps  the  most  obvious  1b  variations  between  the  test 
vehicles  themselves.  The  sharp,  characteristically  super¬ 
sonic  lines  of  the  Basic  Firmer  are  no  doubt  impossible  to 
duplicate  exactly.  Possibly  small  differences  in  construction 
or  even  in  surface  roughness  could  have  great  enough  influ¬ 
ences  on  flow  development  to  affect  materially  the  lift  and 
drag  characteristics  of  the  body  and/or  fins.  A  second  source 
is  errors  in  Betting  and/or  measuring  physical  parameters. 

Even  Bmall  errors  in  such  quantities  as  angle  of  attack  or 
flow  velocity,  or  in  the  alignment  of  balances,  can  produce 
large  discrepancies  in  reported  values  of  force  anu  moment 
coefficients.  A  third  source  common  to  all  experimental 
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techniques.  Is  the  presence  of  the  strut,  l.e.,  the  member 
whloh  supports  the  test  vehicle*  Often,  the  measurements 
made  Include  the  effects  of  forces  and  moments  acting  on  the 
strut  as  well  as  those  acting  on  the  test  vehicle  Itself.  In 
other  cases,  the  measurements  Include  only  the  reactions  ex¬ 
erted  on  the  test  vehicle.  In  either  event,  the  acting  foroea 
and  moments  can  be  different  from  those  which  would  aot  If  the 
strut  were  not  present  In  the  flow  field.  In  certain  of  the 
Basic  Firmer  experiments,2*4  attempts  were  made  to  remove 
strut  effects  from  the  data  by  methods  assuming  linear  super¬ 
position  of  hydrodynamic  effects  of  different  origin.  In  the 
remaining  cases,  strut  interference  either  was  assumed  neg¬ 
ligible,3*9  or  its  influence  was  noted  but  not  evaluated 
quantitatively.9 

The  remaining  factor  contributing  to  the  discrepan¬ 
cies  in  the  hydrodynamic  coefficients  Is  the  presupposition 
of  functional  relationships  between  the  hydrodynamic  forces 
and  moments  and  the  various  motion  parameters.  This  factor 
is  characteristic  of  the  techniques  categorized  here  as 
"oscillatory"  methods  {see  Table  i) ,  These  techniques  depend 
upon  the  assumption  that  the  hydrodynamic  forces  and  moments 
acting  during  the  tests  are  linear  functions  of  the  various 
velocity  and  acceleration  components  of  the  test  vehicle  re¬ 
lative  to  the  fluid.  However,  since  the  test  motions  general¬ 
ly  are  extremely  complicated  oscillatory  motions  (much  more 
complicated  than  the  trajectories  now  being  studied),  it  is 
questionable  whether  the  assumption  is  well  founded  and,  con¬ 
sequently,  whether  the  oscillator  derived  values  for  the 
hydrodynamic  coefficients  are  reliable. 

Table  II  also  Includes  a  row  of  items  labeled, 
"Values  Estimated  From  Theory."  This  row  lists  discrete 
values  for  all  the  required  hydrodynamic  coefficients  except 
the  four  static  and  damping  derivative  coefficients, 

Zw*  Mw*  Zq*  and  Mq*  For  these  quantities,  the  table  refers 


R-898 

-14- 


to  Fig.  3*  whioh  shows  theoretically  derived  ranges  for  their 
values  and  the  interrelationships  between  them. 

The  derivation  of  all  the  results  presented  in 
Table  II  and  Fig.  3  is  described  in  detail  in  Appendix  B, 
However,  the  more  important  features  of  the  derivations  are 
pointed  out  and  disoussed  in  the  present  section. 

The  statio  and  damping  derivatives  were  estimated 
using  two  different  theories.  It  is  inherent  in  both  theories 
that  the  statio  moment  derivative  coefficient,  M^,  and  the 
damping  foroe  and  moment  derivative  coefficients,  Z^  and  M^, 
respectively,  are  all  strongly  dependent  upon  the  statio 
force  derivative  coefficient,  Z^.  In  fact,  the  only  differ¬ 
ences  in  the  results  derived  from  the  two  theories  as  applied 
herein  stem  from  differences  in  the  estimates  for  Z^.  Bach 
method  uses  slender  body  theory  to  predict  the  contribution 
to  Z^  made  by  the  bare  hull.  One  of  the  methods,  due  to 
Flax  and  Lawrenoe,3*  uses  low  aspect  ratio  wing  theory21’  to 
predict  the  contribution  of  the  fins,  and  the  results  of  Low 
and  Stone,84  Spreiter,83  and  Lennarts40  to  account  for  the 
fin-body  interference  effeots.  The  other  method,  used  in  the 
past  by  Davidson  Laboratory  and  others  in  airohip  and  torpedo 
work,  uses  high  aspoct  ratio  wing  theory,  in  conjunction  with 
a  fin  aspect  ratio  augmentation  based  on  geometric  considera¬ 
tions,  to  predict  z{,  the  contribution  to  Z^  due  to  the  force 
on  the  fin  plus  all  fin-body  interactions.  The  low  aspect 
ratio  theory  predicts  a  value  (Z^  -  -13.7)  considerably 
greater  in  magnitude  than  that  derived  from  the  airship  theory 
(Z^  ■  -11.6).  Reference  to  Tabic  II,  however,  reveals  that 
both  values  are  included  in  the  range  of  the  experimental  re¬ 
sults  for  Z^. 

A  review  was  made  of  previous  test  results,  to  com¬ 
pare  experimentally  determined  values  for  zj[  with  values  com¬ 
puted  on  the  basis  of  each  theory.  The  results  of  the  review 
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war*  inconclusive,  as  illustrated  by  the  comparisons  pre¬ 
sented  in  Table  III  for  four  different  hydrodynamlo  mlsailas. 
The  airship  theory  value  is  in  closer  agreement  with  the 
measurement  in  three  of  the  four  oases.  However,  the  opposite 
is  true  in  the  fourth  case,  whioh  is  one  of  a  number  of  such 
oases  reported  by  Plax  and  Lawrence. aa 

The  influence  of  Z^  on  the  other  static  and  damping 
derivatives  is  dependent  upon  the  parameter  xx,  the  longi¬ 
tudinal  coordinate  of  the  point  of  application  of  tx  (see 
Appendix  B) .  No  theoretical  method  for  a  preoise  determina¬ 
tion  of  this  quantity  is  known.  However,  experimental 
evidence1  indicates  that  xj  corresponds  to  a  longitudinal  co¬ 
ordinate  of  the  fins  on  the  body.  Various  values  have  been 
used  for  this  quantity  in  previous  studies.  In  particular, 
values  corresponding  to  the  fin  mid-chord  and  fin  quarter 
S  ohord  (the  aerodynamic  center  in  airfoil  theory)  both  have 

|  been  utilised, 1*,i4  In  the  present  study,  estimates  for  tty 

|  Zq,  and  were  made  on  the  basis  of  xj  values  corresponding 

to  all  longitudinal  coordinates  on  the  fins,  i.e,,  from 
;  -  -0.30,  the  leading  edge,  to  xj  -  -0.40,  the  trailing 

edge.  The  graphs  of  Fig,  3  present  the  results  as  functions 
of  xj.  Two  ourves  of  values  are  shown  for  each  coefficient. 

|  These  correspond  to  results  cased  on  values  for  z[  derived 

;  using  both  the  low  aopect  ratio  theory  and  the  airship  theory. 

|  The  ranges  of  the  theoretically  estimated  values 

|  for  the  static  and  damping  derivative  coefficients  are 

I  -13.7  4  z;  £  -11.6,  -3.83  <  <  -2.02,  -5.48  <  Z^  £  -3.68, 

!  And  -2.19  £  £  -1.19.  These  ranges  are,  respectively,  17#, 

62#,  39#»  And  59#  of  the  nean  calculated  value  for  each  co¬ 
efficient,  In  view  of  the  magnitudes  of  these  ranges,  the 
theoretical  estimates  appear  to  be  of  little  help  in  estab¬ 
lishing  values  for  the  static  and  damping  derivative  coeffi¬ 
cients  for  use  in  trajectory  predictions.  It  will  be  shown 


H-898 

-16- 


later,  however,  that  the  theoretical  relationships  are  valu¬ 
able  Indeed,  because  of  the  fact  that  they  express  In  mathe¬ 
matical  terms  a  basio  physical  truth,  viz.,  that  the  statio 
and  damping  derivatives  are  not  independent  quantities  but 
that  there  exist  definite  relationships  between  them.  The 
significance  of  this  faot  and  the  value  of  its  mathematical 
expression  by  eqs.  B-l  through  B-4  is  shown  in  the  following 
section. 

The  computed  •  stimates  for  the  remaining  hydro¬ 
dynamic  coefficients,  X£,  X^l,  Z£,  Ml,  and  X^,  were  derived 
with  the  guidance  of  theoretical  anchor  empirioal  relation¬ 
ships  published  in  the  literature.  In  all  these  estimates, 
the  effects  of  the  mutual  interference  between  body  and  fin» 
were  assumed  negligible.  The  drag  coefficient,  X£,  was 
estimated  on  the  basis  of  empirical  results  due  primarily  to 
Hoerner.41  The  longitudinal  virtual  mass  coefficient,  X^, 
was  estimated  on  the  basis  of  empirioal  results  presented  by 

Yee-Tak  Yu,8*'*7  The  normal  virtual  mass  coefficient,  Zi, 

w 

and  the  virtual  moment  of  inertia  coefficient,  M^,  were 
estimated  using  equations  derived  on  the  basis  of  geometric 
considerations,  with  the  bare  hull  effects  estimated  using 
the  Lamb  coefficients,18  and  the  results  of  Yee-Tak  Yu8T  used 
again  to  estimate  the  contributions  of  the  fins.  The  second 
order  derivative,  X^,  was  assumed  to  be  equal  to  Z^  on  the 
basis  of  potential  theory.18 

There  are  two  additional  sets  of  hydrodynamic  coef¬ 
ficient  values  presented  in  Table  II.  Each  of  the  "Average 
Experimental  Values"  waB  obtained  by  averaging  the  values  in 
the  rows  above  it,  i.e.,  the  experimentally  derived  values 
discussed  previously.  The  second  set,  labeled,  "NPQ  Analysis-- 
see  section  beginning  on  p,  10 "  was  derived  from  previous 
work  by  Caster, 11  In  that  study,  all  but  one  of  the  hydro- 
dynamic  terms  appearing  in  a  system  of  nonlinear  motion 
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equations  are  epeoified  and  the  remaining  term  (the  damping 
force  derivative  term)  ie  varied  to  find  the  value  giving  the 
beet  overall  agreement  between  three  different  pairs  of  pre¬ 
dicted  and  observed  Basic  Pinner  trajectories.  The  deriva¬ 
tive  coefficients  presented  here  were  evaluated  from  the 
hydrodynamic  input  data  corresponding  to  the  best  fit. 

The  data  presented  in  Table  II  and  Pig.  3  represent 
the  total  body  of  available  information  from  which  the  hydro- 
dynamic  input  data  required  for  the  Basic  Pinner  trajectory 
predictions  could  be  drawn.  With  the  exception  of  the  values 
derived  from  the  results  of  the  NPQ  curve-fitting  analysis, 11 
this  information  was  obtained  without  recourse  to  trajectory 
measurements.  This  is  a  very  important  distinction,  since 
the  practical  value  of  a  trajectory  prediction  technique  re¬ 
quiring  input  data  obtained  from  trajectory  measurements 
surely  is  limited. 

Becaure  of  the  problems  discussed  previously,  how¬ 
ever,  the  selection  of  one  set  of  hydrodynamic  input  data  from 
the  results  obtained  independent  of  trajectory  measurements 
really  could  not  be  made  in  a  logical  manner.  The  available 
experimental  data  were  incomplete  and  imprecise.  The  sources 
of  experimental  error  were  numerous  and  often  impossible  to 
evaluate  quantitatively.  The  values  estimated  by  theoretical 
calculations  also  were  extremely  uncertain.  Consequently,  the 
primary  hope  for  achieving  successful  trajectory  predictions 
was  that  the  computed  trajectories  would  prove  insensitive  to 
variations  in  the  values  used  for  the  hydrodynamic  derivative 
coefficients. 
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TRAJECTORY  PREDICTIONS 

High  Speed  Vertical  Planar  Trajectories 

The  ten  trajectory  predictions  discussed  in  the 
present  section  correspond  to  experiments  in  which  a  model  is 
launched  with  nonzero  horizontal  velocity  and  zero  vertical 
velocity  and  decelerates  throughout  a  trajectory  in  the  verti- 
cal  plane*  The  test  vehicles  range  in  density  from  g  -  0.81 
to  w  -  4.0.  In  eight  of  the  trajectories  (w  >  l),  the  model 

°  **  t  W  v 

falls  under  the  influence  of  gravityj  in  the  other  two  («*  <  1) 

w  ** 

it  rises.  For  a  model  with  g  -  2.0,  the  longitudinal  position 
of  the  CO  is  varied  from  .025  model  lengths  aft  of  the  CB  to 
.025  model  length  forward  of  the  CB.  In  each  run,  the  model 
is  launched  at  high  enough  velocity  so  that  the  flow  is  tur¬ 
bulent  throughout  virtually  the  entire  trajectory.  Hence  the 
high  speed  vertical  planar  trajectories  are  representative  of 
certain  motions  executed  by  full  scale  torpedoes,  submarines, 
and  underwater- launched  missiles. 

The  high  speed  vertical  planar  trajectory  predic¬ 
tions  were  made  using  eqs.  2  through  4.*  The  equations  were 
solved  by  the  numerical  procedure  presented  in  Appendix  C. 

The  initial  conditions  and  other  input  data  correspond  to 
model  trajectory  measurements  supplied  directly  by  the  Bureau 
of  Naval  Weapons,  Since  the  model  is  launched  on  a  horizontal 
tangent,  all  of  the  initial  observed  values  of  velocity  and 
displacement  components  are  zero  except  for  the  initial  longi¬ 
tudinal  velocity,  u(o).  The  value  of  u(o)  is  given  by  the 
slope  of  the  observed  xQ  versus  t  curve,  evaluated  at  t  =  0. 
But  this  slope  cannot  bn  evaluated  precisely  since  the  curve 
is  not  defined  for  t  <  0.  To  circumvent  this  difficulty,  the 
initial  conditions  for  the  trajectory  predictions  were  evalu¬ 
ated  not  at  the  actual  instants  of  launching,  but  at  times 
somewhat  after  launching  where  the  data  curves  are  more  clear¬ 
ly  defined.  Typically,  the  initial  value  of  t  used  in  the 
*  Plus  the  subsidiary  relation,  q  *  9. 
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TRAJECTORY  PREDICTIONS 


High  Speed  Vertical  Planar  Trajectories 

The  ten  trajectory  predictions  discussed  in  the 

present  section  correspond  to  experiments  in  which  a  model  is 

launched  with  nonzero  horizontal  velocity  and  zero  vertical 

velocity  and  decelerates  throughout  a  trajectory  in  the  verti- 

cal  plane.  The  test  vehicles  range  in  density  from  ^  ■  0.81 

to  g  -  4.0.  In  eight  of  the  trajectories  (g  >  l),  the  model 

falls  under  the  influence  of  gravity!  in  the  other  two  (w  <  l) 

w  D 

it  rises.  For  ,  model  wi“h  g  -  2.0,  the  longitudinal  position 
of  the  CO  is  varied  from  .025  model  lengths  aft  of  the  CB  to 
.025  model  length  forward  of  the  CB.  In  each  run,  the  model 
is  launched  at  high  enough  velocity  so  that  the  flow  is  tur¬ 
bulent  throughout  virtually  the  entire  trajectory.  Hence  the 
high  speed  vertical  planar  trajectori'^  are  representative  of 
certain  motions  executed  by  full  scale  torpedoes,  submarines, 
and  underwater-launched  missiles. 

The  high  speed  vertical  planar  trajectory  predic- 

£ 

tions  were  made  using  eqs.  2  through  4.  The  equations  were 
solved  by  the  numerical  procedure  presented  in  Appendix  C. 

The  initial  conditions  and  other  input  data  correspond  to 
model  trajectory  measurements  supplied  directly  by  the  Bureau 
of  Naval  Weapons.  Since  the  model  is  launched  on  a  horizontal 
tangent,  all  of  the  initial  observed  values  of  velocity  and 
displacement  components  are  zero  except  for  the  initial  longi¬ 
tudinal  velocity,  u(o).  The  value  of  u(o)  is  given  by  the 
slope  of  the  observed  xn  versus  t  curve,  evaluated  at  t  ■  0, 
But  this  slope  cannot  be  evaluated  precisely  since  the  curve 
is  not  defined  for  t  <  0.  To  circumvent  this  difficulty,  the 
initial  conditions  for  the  trajectory  predictions  were  evalu¬ 
ated  not  at  the  actual  instants  of  launching,  but  at  times 
somewhat  after  launching  where  the  data  curves  are  more  clear¬ 
ly  defined.  Typically,  the  initial  value  of  t  used  in  the 
*  Plus  the  subsidiary  relation,  q  ■  9. 
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predictions  corresponds  to  a  time  approximately  0.1  sec  after 
launching.  The  values  of  the  initial  conditions  and  all  other 
Input  dat&  for  each  of  the  ten  high  speed  vertical  planar 
trajeotory  predictions  are  given  in  Table  IV. 

The  values  of  the  hydrodynamic  derivative  coeffi¬ 
cients  used  in  the  final  high  speed  vertical  planar  trajeotory 
predictions  were  theoretically  estimated  values  (listed  in 
Table  V) .  In  particular,  the  estimates  for  the  static  and 
damping  derivative  coefficients,  Z^,  M^,  Z^,  and  M^,  are  those 
based  on  the  vai.ue  of  z{  ( zj  «  -9.6)  computed  using  airship 
theory,  and  a  value  of  x|  (x{  -  -0.31)  corresponding  to  a 
polrt  10#  of  the  chord  length  aft  of  the  fin  leading  edge 
(see  Pig.  3).  As  might  be  expected,  the  selection  of  these 
values  for  use  in  the  trajeotory  predictions  was  by  no  means 
a  straightforward  matter.  It  was  made  only  after  a  series  of 
preliminary  calculations  had  been  performed  which  yielded  much 
important  and  surprising  information  regarding  the  influence 
of  the  hydrodynamic  derivative  coefficients  on  the  predicted 
trajectories. 

The  measured  trajectory  chosen  as  the  model  for  tha 
first  preliminary  calculations,  CIT  Run  P-56,  *  is  one  of  the 
trajectories  to  which  prediction  techniques  previously  have 
been  applied  at  NPO  by  Caster.11  In  that  previous  work,  all 
but  one  of  the  hydrodynamic  terms  in  a  system  of  nonlinear 
motion  equations  are  specified  and  the  remaining  term  (the 
damping  force  derivative  term)  is  varied  to  find  the  value 
giving  the  best  agreement  between  computed  and  observed  tra¬ 
jectories,  Although  the  concept  of  varying  this  term  Inde¬ 
pendently  of  the  other  static  and  damping  derivative  coeffi¬ 
cients  1b  not  realistic  since  it  ignores  the  physical  inter- 


*  The  characteristics  of  the  CIT  model  trajectories  are  sum¬ 
marized  in  Table  IV. 
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dependence  of  theee  quantities,*  the  degree  of  acouraoy 
possible  through  using  this  curve- fitting  technique  appears 
significaht.  In  the  three  trajectories  considered  in  the 
previous  study,  the  agreement  between  observed  and  predicted 
values  of  longitude,  depth,  and  inclination,  respectively,  is 
to  within  0.19  model  lengths,  0.12  model  lengths,  and  1.1  dege 
(during  the  course  of  trajectories  of  from  6,6  model  lengths 
to  8.6  model  lengths  in  total  length). 

As  noted  in  the  previous  section,  the  primary  hope 
for  the  suocess  of  the  present  trajeotory  predictions  was  that 
the  computed  trajectories  would  prove  so  insensitive  to 
variations  in  the  hydrodynamic  derivative  coefficients  that 
even  the  wide  ranges  in  the  available  values  for  these  quan¬ 
tities  would  not  introduce  significant  errors  in  the  predic¬ 
tions.  This  hope  was  soon  dispelled  by  calculations  based  on 
two  different  sets  of  values  for  the  hydrodynamic  coefficients 
One  calculation  was  made  using  the  values  evaluated  from  the 
NP3  data  (see  Table  II).  The  Davidson  Laboratory  prediction 
based  on  theee  coefficients  wag  virtually  identical  to  the 
NPO  prediction  and  hence  agreed  well  with  the  observed  re¬ 
sults.  The  other  calculation  waB  made  using  the  "average  ex¬ 
perimental  values"  for  the  coefficients  (see  Table  II) ,  In 
this  case,  the  predicted  trajectory  differed  markedly  from 
the  measured  results.  In  particular,  the  maximum  differences 
between  predicted  and  observed  values  for  depth  and  inclina¬ 
tion  were  about  20  times  as  great  and  about  40  times  as  great, 
respectively,  as  in  the  previous  calculation. 


*  It  is  found  that  the  hydrodynamic  static  and  damping  de¬ 
rivative  coefficients  derived  from  the  NPO  results  do  r.oc 
satisfy  eqs.  B-i  through  J3~4  for  any  reasonable  set  of 
values  for  the  quantities  appearing  therein,  For  this 
reason,  their  physical  validity  is  strongly  doubted. 
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The  results  of  these  two  calculations  demonstrate 
conclusively  that  the  scatter  In  the  available  values  for  the 
hydrodynamic  derivative  coefficients  was  significant  from  the 
point  of  view  of  its  effect  on  predicted  trajectories.  Ad¬ 
ditional  preliminary  calculations  substantiated  this  result, 
and  also  yielded  information  on  the  basis  of  which  further 
conclusions  could  be  drawn.  In  all,  a  great  many  preliminary 
calculations  were  made,  and  their  results,  taker  individually, 
have  little  meaning.  Accordingly,  these  resultB  will  not  be 
presented  in  detail.  Instead,  the  major  conclusions  to  which 
they  led  will  be  discussed,  and  the  particular  resultB  lead¬ 
ing  to  these  conclusions  will  be  noted. 

It  was  found  that  the  predicted  trajectories  were 
extremely  sensitive  to  independent  variations  of  the  hydro- 
dynamic  static  and  damping  derivative  coefficients.  For  ex¬ 
ample,  a  mere  10#  change  in  the  value  of  --  much  less  than 
the  range  of  available  values  for  this  coefficient  (Bee 
Table  II)  was  found  to  produce  changes  in  predicted  values  of 
depth  and  inclination  as  great  sb  0.19  model  lengths  and  2.9 
degs,  respectively.  These  variations  are  both  considerably 
greater  than  the  maximum  differences  between  observed  and  com¬ 
puted  results  obtained  in  the  NPG  study.11  As  noted  previous¬ 
ly,  however,  independent  variations  of  the  static  and  damping 
derivative  coefficients  are  not  realistic  in  a  physical  sense. 
Accordingly,  the  effect  on  the  predicted  trajectories  of  re¬ 
stricting  the  static  and  damping  coefficients  to  mutually 
consistent  combinations  was  investigated,  by  comparing  the 
results  of  trajectory  predictions  based  on  theoretically 
estimated  values  for  these  coefficients  (see  Fig.  3)  with 
those  of  predictions  based  on  coefficients  obtained  from  other 
sources.  Since  all  theoretical  values  were  derived  by  using 
eqB.  B-l  through  B-4,  they  are  necessarily  self-consistent. 

The  coefficient  values  obtained  f’rom  other  sources  generally 
are  not It  was  found  that  the  sensitivity  of  the  computed  tra¬ 
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Jectories  to  variations  in  the  values  used  for  the  derivatives 
was  reduced  drastically  by  restricting  the  combinations  to 
satisfy  the  self-consistency  condition.  The  graphs  in  the 
left-hand  column  of  Pig.  4  show  the  maximum  differences 
between  observed  and  computed  valuos  of  longitude,  depth,  and 
Inclination  obtained  in  predictions  for  CIT  Run  F-56  made 
using  theoretical  estimates  for  the  static  and  damping  de¬ 
rivative  coefficients  computed  on  the  basis  o r  x{  values  from 
-0,30  to  -0.35  and  zi  values  of  -9.6  and  -11.7  (see  Pig.  3). 

The  values  for  the  other  hydrodynamic  coefficients  used  in 
all  these  calculations  are  the  theoretical  estimates  given  in 
Table  II.  It  is  found  that  the  variations  in  the  plotted 
values  from  one  calculation  to  another  are  generally  small, 
particularly  compared  with  differences  attained  in  many  of  the 
calculations  in  which  inconsistent  values  were  uBed  for  the  static 
and  damping  derivative  coefficients.  An  example  of  the  re¬ 
sults  of  this  latter  type  is  represented  on  the  graphs  by  the 
values  obtained  on  the  basis  of  the  average  measured  coeffi¬ 
cients. 

It  also  is  seen  from  the  graphs  of  Pig.  4  that  dif¬ 
ferences  between  the  observed  results  and  the  Davidson  Labora¬ 
tory  preliminary  calculations  for  CIT'  Run  P-56  are  generally 
of  the  same  order  of  magnitude  as  those  obtained  in  the  pre¬ 
vious  NPG  curve-fitting  analysis.  This  result  gave  rise  to 
new  optimism  regarding  the  prediction  of  the  high  speed 
vertical  planar  trajectories,  and  prompted  extending  the  pre¬ 
liminary  calculations  to  additional  model  trajectories.  The 
results  of  these  calculations  are  summarized  by  the  graphs  in 
the  rl^ht-hand  column  of  Fig.  4,  which  present  data  similar  to 
that  in  the  left-hand  column  graphs  discussed  previously,  but 
on  the  bar  is  of  average  values  for  three  runs,  CIT  Nos.  F-56, 
F-59,  anu  P-64,  The  average  magnitudes  of  the  differences 
between  observed  and  predicted  values  generally  are  even 
smaller  than  the  corresponding  quantities  based  on  the  P-56 
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prediction  alone.  In  addition,  the  variations  in  the  average 
quantities  with  x{  and  z{  are  smaller.  On  the  basis  of  tneoe 
results,  it  was  concluded  that  reasonably  accurate  trajectory 
predictions  oould  be  obtained  on  the  basis  of  theoretically 
estimated  values  for  the  static  and  damping  derivative  coef¬ 
ficients  computed  using  any  x{  and  z{  values  in  the  ranges  in¬ 
cluded  in  Pig,  4.  On  the  average,  however,  the  closest  agree¬ 
ment  with  experiment  seemed  to  be  achieved  using  the  z[  value 
oomputed  using  airship  theory  (z{  -  -9.6)  and  the  value 
xl  -  -0.31.*  It  is  for  this  reason  that  these  values  were 
selected  to  be  used  in  the  subsequent  trajectory  predictions. 

The  selection  of  the  theoretically  estimated  values 
for  the  remaining  required  coefficients  was  primarily  in  the 
Interest  of  consistency.  In  any  event,  in  the  case  of  the 
drag  coefficient,  X^,  the  theoretical  value  is  identical  to 
the  average  measured  value  and  hence  1b  the  only  logical 
choice.  In  the  case  of  the  remaining  coefficients  (the 
virtual  mass  and  moment  of  inertia  coefficients),  the  question 
is  an  academic  one  since  the  results  of  the  preliminary  calcu¬ 
lations  indicated  that  any  reasonable  variations  in  these 
values  have  but  negligible  effects  on  the  predicted  trajec¬ 
tories.** 


*  The  selection  of  this  value  for  x]  was  also  influenced  by 
the  results  of  preliminary  calculations  corresponding  to 
CIT  Run  P-63  which  were  not  included  in  Pig.  4,  because 
they  did  not  cover  the  complete  range  of  x{  and  z{  values 
considered  there.  However,  the  P-63  results  do  not  differ 
qualitatively  from  those  for  CIT  Runs  P-56,  F-59,  and  F-64 
discusaed  above. 

**  This  result  is  consistent  with  the  results  of  an  analytic 
variation  of  parameter  analysis  made  later  for  the  case  of 
the  straight  horizontal  trajectory. 


Curves  of  the  values  computed  for  xQ,  zQ,  xQ,  z0, 
and  0  In  each  of  the  final  high  speed  vertical  trajectory 
predictions  are  plotted  versus  t  In  Figs.  5  through  14, 
together  with  the  corresponding  values  observed  In  the  model 
experiments.  Tabulations  of  the  extreme  values  of  the  dif¬ 
ferences  between  the  observed  and  computed  values  of  each  of 
these  quantities  are  presented  in  Table  VI.  The  predicted 
values  agree  with  the  observations  (during  the  course  of 
trajectories  of  from  6.6  to  9.6  model  lengths  in  total  length) 
to  within  0.22  model  lengths  in  longitude,  0.07  model  lengths 
in  depth,  and  1.8  dogs  in  inclination.  In  terms  of  percent¬ 
ages  of  maximum  measured  values,  these  figures  correspond  to 
less  than  3$  in  longitude,  4$  in  depth,  and  6%  in  inclination, 
which  indicates  that  the  predictions  are  successful  from  an 
engineering  point  of  view.  Relative  to  the  results  obtained 
in  the  earlier  NPO  curve-fitting  analysis,  the  present  re¬ 
sults  also  are  impressive.  Ab  can  be  seen  from  Table  VII. 
which  compares  the  Davidson  Laboratory  and  NPO  results  for 
all  three  trajectories  considered  in  the  earlier  work,  the 
extreme  values  of  the  differences  between  observed  and  com¬ 
puted  trajectory  parameters  in  most  cases  are  smaller  for  the 
present  predictions  than  for  the  corresponding  NPO  results. 

A  significant  feature  of  the  results  of  the  high 
speed  vertical  planar  trajectory  predictions  is  the  relative 
agreement  with  experiment  achieved  by  those  predictions  in¬ 
volved  in  the  preliminary  calculations  and  those  which  were 
not.  As  noted  previously,  the  values  of  the  hydrodynamic 
static  and  damping  derivative  coefficients  used  in  the  final 
predictions  were  e  itablished  with  the  guidance  of  the  results 
of  CIT  Runs  F-56,  F-59>  F-63»  and  F-64  (Davidson  Laboratory 
calculations  2,  5,  9*  and  10) .  The  remainder  of  the  high 
speed  vertical  planar  runs  were  not  involved  in  this  work. 

As  can  be  seen  from  Table  VI  and  Figs,  5  through  14,  however, 
the  agreement  with  experiment  achieved  by  a  prediction  is  not 
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noticeably  dependent  upon  whether  or  not  it  was  involved  in 
the  preliminary  calculations.  This  is  regarded  as  strong 
verification  of  the  validity  of  the  prediction  techniques 
which  were  employed. 

However,  the  opposite  result  was  found  for  predic¬ 
tions  made  using  the  hydrodynamic  coefficients  evaluated  from 
the  NPO  data.  The  preliminary  calculations  included  pre¬ 
dictions  made  using  these  values,  for  two  trajectories  (CIT 
Runs  P-59  and  P-63)  which  had  not  been  considered  previously 
in  the  NPO  analysis.  In  both  cases,  differences  between  ob¬ 
served  and  predicted  values  01  both  depth  and  inclination 
angle  were  obtained  which  were  greater  ( often  more  than  twice 
as  great)  than  the  maximum  differences  obtained  for  all  three 
runs  considered  in  the  earlier  study.  This  result  was  not 
surprising,  since  the  failure  of  the  hydrodynamic  static  and 
damping  derivative  coefficients  obtained  from  the  NPO  data  to 
satisfy  the  self-consistency  condition  implied  that  they  were 
not  physically  realistic  and,  consequently,  that  predictions 
based  on  them  should  not  be  expected  to  be  generally  success¬ 
ful. 


In  addition  to  the  numerical  solution,  which  mav  be 
regarded  as  exact,  an  analytic  approximate  solution  applic¬ 
able  to  the  high  speed  vertical  planar  trajectories. was  de¬ 
veloped.  The  approximate  oolution  is  valid  for  trajectories 
in  which  the  transverse  motions  are  not  excessive.  In  parti¬ 
cular,  it  is  assumed  that  the  following  conditions  are  satis¬ 
fied  throughout  the  trajectory: 


(W-B)  sine 

!  W* 


w2 


q2  -  wq  s  0 


cos  0=1 


(9) 

(10) 


The  approximate  solution  obtained  on  the  basis  of 
these  assumptions  is  presented  in  Appendix  C.  The  results  of 
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the  approximate  solutions  for  four  typical  cases  and  the  cor¬ 
responding  results  of  the  numerical  solutions  are  presented 
in  Fig.  15.  The  approximate  and  exact  solutions  are  in  good 
agreement  except  near  the  end  of  each  trajectory.  Therefore, 
the  approximate  solution  is  of  value  on  two  counts.  First  of 
all,  it  enables  the  engineer  to  make  reasonably  good  trajec¬ 
tory  predictions  with  relatively  little  work.  Secondly,  it 
lends  itself  to  straightforward  interpretation  and  a  direot 
method  of  estimating  the  influence  on  the  motion  of  the  var¬ 
ious  parameters  which  appear  in  the  motion  equations.  In  the 
latter  connection,  it  permitB  the  problem  to  be  treated  in  a 
manner  which  is  not  possible  when  only  a  numerical  solution 
is  available.  However,  since  the  variation  of  parameter 
analysis  is  extremely  complicated,  it  has  not  been  Included 
in  the  present  study.  The  general  method  is  illustrated, 
however,  by  the  analysis  performed  later  for  the  much  simpler 
case  of  the  straight  horizontal  trajectory. 


Straight  Horizontal  Trajectory 


The  trajectory  prediction  discussed  in  the  present 
section  corresponds  to  an  experiment  in  which  a  neutrally 
buoyant  model,  whose  center  of  gravity  and  center  of  buoyancy 
coincide,  is  launchod  in  a  horizontal  direction  and  deceler¬ 
ates  throughout  a  one-dimenBional  trajectory  along  the  line 
of  launching.  The  trajectory  actually  is  a  special  case  of 
the  high  speed  vertical  planar  motions  discussed  earlier  in 
the  section  beginning  on  p.  19.  It  is  considered  here  separ¬ 
ately,  however,  because  restriction  to  straight  horizontal 
motion  permitB  use  of  eq.  6,  which  may  be  integrated  directly 
to  yield  the  solution 


u 


u(0) 


“  bf"+"  1 


(11) 
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where, 

fw<°> 

"  ■ 

Since  i.Q  -  u  for  the  straight  horizontal  trajectory, 

j.  _  u(  0) 

*o  “  BT  +  T 

The  solution  for  the  longitudinal  position  ie  ob¬ 
tained  by  integrating  eq,  13  and  making  use  of  the  initial 
condition  to  get 

x0  -  log  (bfc  +  l)  +  x0(°)  (l1*) 

The  initial  conditions  and  other  input  parameters 
for  the  straight  horizontal  trajectory  prediction  are  given 
in  Table  IV.  Again,  the  initial  value  of  time  used  in  the 
calculation  does  not  correspond  to  the  instant  of  launching, 
but-  rather  to  a  time  0.02  sec  after  the  launching  ( see  dis¬ 
cussion  beginning  on  p.  19  ). 

Figure  16  shows  a  comparison  of  the  observed  and 
predicted  values  of  velocity  and  longitude  as  functions  of 
time.  As  is  the  case  for  the  more  general  high  speed  vertical 
trajectory  predictions,  the  computed  values  are  in  substantial 
agreement  with  the  experimental  data.  Quantitatively,  the 
agreement  in  longitude  is  to  within  0.35  model  lengths,  or 
less  than  3.1$  of  the  total  trajectory  length. 

Since  an  analytic  solution  is  available,  it  is 
possible  to  make  a  straightforward  analysis  to  assess  the 
effects  on  the  predicted  trajectories  of  possible  errors  in 
the  values  of  the  input  data.  Foi  simplicity,  the  present 


(12) 

(13) 
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analysis  considers  the  effects  of  errors  in  only  the  most  ob¬ 
vious  sources,  the  hydrodynamio  coefficients  and  the  initial 
conditions.  On  this  basis,  the  total  differential  of  the 
longitudinal  velocity  is  given  by 

S*  a* 

"  SJTJ  "o  +  ^mT  dm*  +  SOT  du(0)  '  (15) 

where  the  apparent  mass, 

mx-  m  -  §AlXt  (16) 

C  A  U 


has  been  introduced  to  add  symmetry  to  the  resulting  express¬ 
ions.  Substituting  eq.  13  into  eq.  15  and  evaluating  the  de¬ 
rivatives  that  appear  in  the  resulting  equation  gives 


X_  ..  x„ 

o  jv  i  bt  o  i  _  o. 


'  stn  n  dx«  +  bht  %  'Jm>  +ut§)(1'  wrr)du(o) 

0  (17) 


whence,  to  vhe  first  order  of  approxjmation. 


It  may  be  concluded  that  inaccuracies  in  the  pre¬ 


dicted  results  due  to  errors  in  evaluating  u(o)  are  not  sig¬ 
nificant.  The  basis  for  this  conclusion  is  that  the  time- 
dependent  error  term  which  is  proportional  to  Au(o)  never 
becomes  aB  large  as  the  constant  term,  which  is  known 

to  be  insignificant  because  of  the  perfect  agreement  between 
the  observed  and  predicted  values  at  the  beginning  of  the 
trajectory. 


To  evaluate  the  relative  sensitivity  to  errors  in 

and  X^J,  data  from  Tables  IV  and  V  are  substituted  into 

eq.  16  and  the  result  is  compared  with  the  value  of  X!.  to 

u 

find 
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mi«.  X^J  -  -0.454x^  (19) 

But,  on  the  basis  of  the  assumptions  of  the  analysis, 

Am1— |AxMXl  (20) 

whence,  substituting  from  Table  IV, 

Amx  -  -  (.0211)  (1.667)  AX!  .0352  AX^  (2l) 

Substituting  eqs.  19  and  21  into  eq.  18,  and  assuming  on  the 
basis  of  the  preceding  argument  that  the  errors  proportional 
to  Au(o)  are  negligible,  gives 


A*o 


(AX!  AX’  \ 

,0775  T"^/ 


bt 

FETT 


(22) 


Prom  this  result.  It  Is  seen  that  an  error  of  a 


certain  percentage  In  X£  has  an  effect  on  the  predicted  value 
of  almost  thirteen  times  as  great  as  does  an  error  of  the 
aims  percentage  In  Xh 

Further,  it  Is  noted  that,  for  t  >  2  sec,  the 
approximation 


bt  s 
bt+1 


(23) 


is  accurate  to  within  better  than  10#.  Substituting  this  re¬ 
sult  into  eq.  22  gives 


AX!  AX’ 
3  (.0775  -  -j?) 

U  c> 


(24) 


Interpreting  AxQ  to  be  the  difference  between  the 
observed  and  predicted  valuee  of  longitudinal  velocity,  the 
trajectory  data  for  t  >  2  sec  are  used  In  conjunction  with 
eq.  24  to  give,  approximately, 
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where  X'  +  AX'  and  XI  +  AXI  are  Identified  aa  values  of  the 
oo  u  u 

drag  coefficient  and  the  added  mass  coefficient,  respectively, 
which  would  yield  predictions  in  perfeot  agreement  with  the 
observed  values. 


Obviously,  there  are  an  infinite  number  of  oombina- 
AX!  AX' 

tions  of  -ft  and  *^r  satisfying  eq.  25.  The  most  interesting 

of  these  are  the  cases  where  all  of  the  error  is  attributed 
to  either  one  or  the  other  of  the  coefficients.  In  particular, 
it  is  indicated  that  either  a  10£  increase  in  the  magnitude 
of  X^  or  a  130 f  decrease  in  the  magnitude  of  X^  would  bring 
the  observed  and  predicted  results  into  agreement.  The  second 
alternative  may  be  rejected  as  physically  unreasonable  ainoe 
it  oorrespords  to  a  negative  value  for  the  added  mass.  The 
first  alternative,  however,  is  reasonable  when  considering 
the  available  experimental  data  for  X£  (see  Fig.  2). 


Low  Speed  Vertical  Planar  Trajectories 

The  eight  trajectory  predictions  discussed  in  the 
present  section  correspond  to  experiments  in  which  a  heavier- 
than-water  model  is  dropped  from  rest  with  an  initial  value 
of  pitch  angle  not  equal  to  -90°,  and  accelerates  throughout 
a  trajectory  in  the  vertical  plane.  The  test  vehicles  vary 
in  density  from  W/B  -  1.1  to  W/B  *  4.0.  For  a  model  with 
W/B  ■  2.0,  the  longitudinal  position  of  the  CQ  is  varied  from 
.005  model  lengths  aft  of  the  CB  to  .012  model  lengths  forward 
of  the  CB.  Since  the  models  are  dropped  from  rest,  the 
motions  run  in  or,  at  least,  through  the  laminar  flow  regime. 
Therefore,  since  actual  motions  of  full-scale  missiles  are 
invariably  in  the  turbulent  flow  regime,  the  low  speed 
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vertical  planar  trajectory  predictions  are  of  little  practi¬ 
cal  interest.  In  fact,  their  results  are  reported  largely 
because  of  their  remarkable  accuracy  even  though  they  were 
derived  on  the  basis  of  methods  chat  seeded  grossly  inappli¬ 
cable.  The  first  of  these  was  the  use  of  motion  equations 
(eqs.  2  through  4),  formulated  for  conditions  which  are  not 
satisfied  throughout  by  the  low  speed  vertical  planar 
trajectories.  This  question  has  already  been  discussed  at 
length,  in  the  second  section.  The  second  was  the  use  of  the 
same  constant  values  for  the  hydrodynamic  coefficients  as 
were  employed  in  the  high  speed  trajectory  predictions  (see 
Table  V) .  This  step,  which  is  contrary  to  accepted  theory,1® 
was  necessitated  by  the  limitations  of  the  available  hydro- 
dynamic  data  ( see  Table  I)  and  the  lack  of  reliable  techniques 
for  theoretical  estimation  of  the  behavior  of  the  hydrodynamic 
derivative  coefficients  at  low  Reynolds  numbers. 

The  initial  conditions  and  other  input  parameters 
for  the  low  speed  vertical  trajectory  predictions  correspond 
to  model  trajectory  tests  reported  by  Price.7  The  initial 
values  of  time  used  in  the  calculations  correspond  to  the  in¬ 
stants  of  dropping  in  the  tests.  Hence  all  of  the  initial 
values  except  0(o)  are  equal  to  zero.  All  input  data  for 
each  of  the  low  speed  vertical  trajectory  predictions  are 
given  in  Table  IV,  Because  of  the  lack  of  practical  interest 
In  the  low  speed  vertical  planar  trajectories,  detailed  re¬ 
sults  of  each  one  of  the  calculations  are  not  presented.  For 
illustrative  purposes,  curves  of  the  values  predicted  for 
x  ,  z  ,  x  ,  z  ,  and  9  in  three  typical  calculations  are 
plotted  versus  t  in  Figs.  17  through  19,  together  with  the 
corresponding  values  observed  in  the  model  tests.  Extreme 
values  of  the  differences  between  the  observed  and  computed 
values  of  xQ,  zQ  and  9  for  each  trajectory  are  presented  in 
Table  VI. 
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The  low  speed  vertical  planar  trajectory  predic¬ 
tions  are  considerably  less  accurate  than  are  the  high  speed 
predictions.  The  average  magnitudes  of  the  maximum  differ¬ 
ences  between  computed  and  observed  values  of  longitude, 
depth,  and  inclination  are,  respectively,  0.14  model  lengths, 
0.14  model  lengths,  and  2.9  dega  (see  Table  VI) .  The  cor¬ 
responding  figures  for  the  high  speed  vertical  trajectory 
predictions  are  0.09  model  lengths,  0.03  model  lengths,  and 
0.9  degs,  respectively  (the  comparison  is  all  the  more  sig¬ 
nificant  since  the  high  speed  trajectories  have  an  average 
total  length,  in  model  lengths  of  more  than  twice  that  of  the 
low  speed  runs) .  However,  the  extreme  value  of  the  differ¬ 
ences  between  the  predicted  and  observed  values  of  depth  has 
a  magnitude  of  but  0.35  model  lengths.  Roughly,  this  means 
that  the  displacement  is  predicted  to  within  less  than  10# 
of  the  maximum  measured  depth.  This  constitutes  an  adequate 
estimate  for  many  engineering  applications,  and  a  remarkably 
accurate  one  in  view  of  the  assumptions  made  in  the  analysis. 

Straight  Vertical  Trajectories 

The  seven  trajectory  predictions  discussed  ih  this 
section  correspond  to  experiments  in  which  a  heavier-than- 
water  model  is  dropped  from  rest,  with  an  initial  pitch  angle 
equal  to  -90°,  and  accelerates  throughout  a  one-dimensional 
vertical  trajectory.  The  trajectories  are  special  cases  of 
the  low  speed  vertical  planar  trajectories  discussed  earlier, 
in  the  section  beginning  on  p.  31.  They  are  considered  here 
separately,  however,  because  of  their  relative  simplicity 
(see  Qiscussion  in  the  first  and  second  sections).  The  pre¬ 
dictions  were  made  using  eq.  8,  which  may  be  written  in  the 
form 

u  -  C2  -  rf  (26) 
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m  -  fAx«I 


na 


m  -  fAx«  J 


(28) 


Note  that  the  aubBtltutlons  are  made  In  such  a  way  that  the 
constants  C  and  rj  are  real  numbers.  Equation  26  is  Inte¬ 
grated  to  obtain 

u  ■  ^  tanh  ( CT)t)  (29) 

where  the  initial  condition  u(0)  -  0  is  incorporated  in  the 
solution.  Therefore,  using  the  Table  of  Direction  Cosines 
(Appendix  A): 

£0  -  i  tanh  (  CT.t)  ( 30) 

The  solution  for  the  depth  is  obtained  by  integrating  eq>  30 
and  using  the  initial  condition,  zQ(o)  -  0  ,  to  get 

zQ  *  log  cosh  (  Ctj t )  (31) 

• 

The  input  data  for  the  straight  vertical  drop  tra¬ 
jectories  are  design  values  (i.e.,  they  correspond  to  the 
standard  Basic  Pinner  configuration)  rather  than  actual 
values  measured  for  the  various  test  vehicles.  The  reason 
for  this  is  that  the  actual  values  of  the  input  parameters 
are  reported  in  ref.  7  for  only  some  of  the  models  tested, 
and  it  was  desired  that  the  predictions  be  made  in  a  consis¬ 
tent  manner.  Preliminary  calculations  showed,  however,  that 
the  magnitudes  of  the  differences  between  the  design  values 
and  those  actual  values  which  are  known  are  not  great  enough 
to  affect  the  predicted  motions  appreciably  except  in  cases 
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where  the  specific  gravity  Is  close  to  unity.  The  values 

used  for  the  two  required  hydrodynamic  coefficients#  and 

Xi#  are  the  constant  theoretical  estimates  derived  as  shown 
u 

In  Appendix  B.  The  values  of  all  of  the  Input  parameters 
for  each  of  the  seven  straight  vertical  trajectory  predic¬ 
tions  are  given  in  Table  IV, 

Curves  of  the  values  computed  for  zQ  and  zQ  in  each 
of  the  straight  vertical  trajectory  predictions  are  plotted 
versus  t  in  Pigs.  20  through  22#  together  with  the  corres¬ 
ponding  values  observed  in  the  model  testa.  It  is  seen  that 
the  predicted  and  observed  values  agree  well  in  every  case. 
With  the  exception  of  Calculation  20,  the  extreme  value  of 
the  difference  between  the  predicted  and  observed  values  of 
depth  has  a  magnitude  of  0.13  model  lengths#  or  about  1.6# 
of  the  maximum  measured  depth  (see  Table  VI) .  The  value  tor 
Calculation  20  is  0.38  model  lengths#  or  about  4.7#  of  the 
total  depth.  But  this  calculation  corresponds  to  experiments 
in  which  certain  modifications,  too  slight  to  be  reflected 
in  the  estimates  of  the  hydrodynamic  force  and  moment  coef¬ 
ficients  were  made  to  the  teBt  vehicles.  Moreover#  the 
models  in  two  of  these  runs  were  released  by  hand  rather  than 
by  the  mechanical  launcher  used  in  the  other  tests.  Hence 
the  greater  discrepancy  was  expected. 

In  spite  of  the  excellent  agreement  with  observed 
motions  achieved  by  the  straight  vertical  trajectory  predic¬ 
tions  made  on  the  basis  of  a  constant  drag  coefficient,  it 
was  felt  that  some  attempt  should  be  made  to  reconcile  the 
results  with  accepted  theory#  i.e.,  to  show  that  the  inclu¬ 
sion  of  a  realistic  Reynolds  number  dependent  component  of 
the  drag  coefficient  does  not  significantly  affect  the  pre¬ 
dicted  trajectories.  Accordingly#  trajectory  predictions 
corresponding  to  CIT  vertical  drops  F-9  and  F-47,  48 
(Davidson  Laboratory  Calculations  No.  21  and  26)  were 
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repeated  on  the  basis  of  Reynolds  number  dependent  drag  coef¬ 
ficients.  These  runs  represent  the  extreme  values  (least  and 
greatest,  respectively)  of  Reynolds  number  variation  In  all 
the  trajectories  of  this  type.  Hence,  It  is  reasoned  that 
any  significant  Reynolds  number  dependent  phenomenon  must 
manifest  itself  in  at  least  one  of  these  two  runs. 


Reynolds  numbor  dependent  drag  coefficients  were 
derived  by  amalgamating  available  Basic  Firmer  data  and  em¬ 
pirical  and  theoretical  results  for  the  low  Reynolds  number 
range  which  are  presented  by  Hoerner.41  Two  different  ex¬ 
pressions  were  employed,  one  based  on  the  laminar  flow  drag 
characteristics  of  a  sphere  and  the  other  on  the  characteris¬ 
tics  of  a  flat  plate.  The  expressions  employed  were  tailored 
to  approximate  the  laminar  flow  data  in  the  low  Reynolds  num¬ 
ber  range,  and  to  approach  the  value  used  in  the  constant 
coefficient  predictions  as  the  Reynolds  number  increases  into 
the  turbulent  flow  regime. 

The  expression  for  based  on  the  characteristics 
of  a  sphere  is 


(32) 


where  Dfl  -  0.45  and  Rfl  -  24.  £  ft/sec.  Figure  23  includes  a 

graph  showing  the  variation  of  X'  with  Reynolds  number  over 

°s 

most  of  the  range  covered  by  the  vertical  drop  trajectories. 

When  eq.  32  for  X'  is  substituted  in  the  straight 
°s 

vertical  trajectory  motion  equation  (eq.  26),  the  result  is 
written  in  the  form 

«  -  u2  -  $2u  (33) 
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5Ax  Da 

(34) 

» -  fv*d 

(35) 

"  -  &c«d 


and  C*  1>  ae  defined  by  eq.  27.  Note  that  the  eubetltutlona 
are  made  In  auoh  a  way  that  the  oonatanta  are  real  numbers. 
Equation  33  la  Integrated  to  obtain 


where 


and  where  the  initial  condition  u(o)  -  0  le  Incorporated  in 
the  solution  and  eq.  C-9  is  used  with  9  ■  -  jy.  The  solution 
for  the  depth  is  obtained  by  integrating  eq.  36  and  using  the 
initial  condition  zQ(o)  -  0,  to  get 


The  expression  for  based  on  the  characteristics 
of  a  flat  plate  is 


a 


Here, 


Rp  -  59.5 


[i  *Vsec] 


V 


but  Dp  is  given  dif¬ 


ferent  values  for  the  two  different  trajectory  predictions. 


In  the  prediction  corresponding  to  CIT  Run  No.  F-9,  D  -  0.35. 


As  seen  from  Fig.  23,  this  value  is  such  that  X’  approaches 

P 


the  constant  coefficient  value  of  -0.45  as  the  Reynolds  num¬ 
ber  approaches  the  maximum  value  attained  in  the  observed 
trajectory  (approximately  0.4  x  10e) . 


In  the  prediction  corresponding  to  CIT  Runs  F-47, 

48,  Dp  -  0.45.  As  seen  from  Fig.  23,  this  value  is  such  that 

X'  is  considerably  greater  in  magnitude  than  the  constant 
P 

coefficient  value  of  -0.45  throughout  the  Reynolds  number 
range  attained  in  the  observed  trajectory.  Nevertheless,  it 
is  sufficient  to  show  that  including  a  term  of  this  type  is 
unimportant,  since  it  overestimates  rather  than  underestimates 
the  effeot. 


No  analytical  solution  was  obtained  for  the  straight 

vertical  trajectory  motion  equation  with  X*  ■  X'  .  It  was 

P 

solved  by  a  numerical  procedure  similar  to  that  used  in  the 
solution  of  the  vertical  planar  trajectory  motion  equations. 

The?  straight  vertical  trajectory  predictions  ob¬ 
tained  on  the  basis  of  the  Reynolds  number  dependent  drag 
coefficients  substantiate  the  predictions  obtained  assuming 
that  the  coefficient  is  constant.  In  none  of  the  four  repeat 
predictions  were  the  results  appreciably  different  from  the 
original  computations#  This  result  indicates  that  terms  pro¬ 
portional  to  accelerations  (primarily  inertial  terms)  are  so 
dominant  in  the  low  velocity  portions  of  the  trajectories 
that  even  large  errors  in  the  hydrodynamic  coefficients  are 
unimportant. 


CONCLUSIONS 


Observed  trajectories  could  not  be  predicted  using 
hydrodynamic  coefficients  determined  from  available  experi¬ 
mental  data  alone.  However,  coefficients  obtained  on  the 
basis  of  existing  theory  yielded  predictions  that  agreed  well 
with  the  observed  motions  (see  Table  VI) . 

Computed  trajectories  were  particularly  sensitive 
to  independent  ohanges  in  the  values  used  for  the  hydrodynamic 
statio  and  damping  derivative  coefficients,  Z^,  M^,  Z^,  and 
M^.  However,  when  the  values  were  kept  consistent  with  re¬ 
lations  expressing  the  physical  interdependence  of  the  coef¬ 
ficients  (eqs.  B-l  through  B-4),  they  could  be  varied  over 
fairly  wide  ranges  without  appreciably  affecting  the  predic¬ 
tions.  This  result  is  particularly  significant,  and  efforts 
should  be  made  to  generalize  it.  In  particular,  its  validity 
for  the  oase  of  a  more  representative  hydrodynamic  missile 
(c.g.,  a  standard  torpedo  or  submarine  type  configuration) 
should  be  investigated.  Such  an  investigation  could  be  in 
the  form  of  calculations  similar  to  those  performed  in  the 
present  study.  The  problem  also  could  bs  approached  analy¬ 
tically,  however,  through  the  use  of  approximate  methods 
which  make  analytic  solutii  n  of  the  motion  equations  possible 
(e.g.,  the  approximate  solution  developed  herein). 

It  is  l'elt  that  the  present  results  demonstrate 
clearly  the  potential  value  of  simple  analytic  methods  for 
predicting  hydrodynamic  derivative  coefficients,  and  it  is 
urged  that  greater  emphasis  be  placed  on  their  use  in  the 
future. 

Prom  a  practical  standpoint,  the  important  results 
of  the  present  study  are  those  pertaining  to  high  speed  tra¬ 
jectories,  The  low  speed  trajectory  predictions  also  are  in¬ 
teresting,  however,  because  of  their  surprising  accuracy. 
Predict  Iona  of  straight  (zero  attack  angle)  vertical  drops 
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from  rest  agreed  well  with  observed  trajectories,  even  though 
a  oonetant  drag  coefficient  corresponding  to  turbulent  flow 
was  used  in  the  oaloulationa.  Low  speed  vertical  planar  tra¬ 
jectory  predictions,  although  markedly  less  accurate  than  the 
high  speed  predictions,  nevertheless  yielded  results  con¬ 
sidered  adequate  for  many  engineering  applications,  in  spite 
of  the  fact  that  the  calculations  took  into  aooount  neither 
nonlinear  hydrodynamio  normal  forces  and  pitching  moments 
(although  angles  of  attack  as  great  as  30°  were  attained)  nor 
variations  of  hydrodynamic  coefficients  with  Reynolds  number. 
The  reason  for  the  unrealistic  simplicity  of  the  mathematical 
model  used  in  the  low  speed  planar  predictions  was  the  laok 
of  any  reliable  experimental  or  theoretical  basis  for  estima¬ 
ting  the  extensive  hydrodynamio  input  data  required  for  a 
more  sophisticated  formulation.  The  surprising  aoouracy  of 
the  predictions  is  attributed  to  the  dominance  of  terms  pro¬ 
portional  to  aooeleration  (i.e.,  inertia!  terms)  in  the  low 
velocity  portions  of  the  trajectories.  It  is  felt  that  any 
improvement  in  the  acouraoy  of  the  low  speed  vertical  planar 
trajeotory  predictions  requires  the  performance  of  experi¬ 
ments  to  obtain  detailed  hydrodynamio  force  and  moment  data 
for  the  Basio  Pinner  in  the  low  Reynolds  number  high  attack 
angle  regime. 

All  available  Basic  Pinner  trajectory  test  results 
have  not  been  considered  in  the  present  study.  In  particular, 
several  heretofore  unpublished  six  degs  of  freedom  trajec¬ 
tories  also  are  available. ae  It  is  recommended  that  attempts 
be  made  to  prediot  these  more  general  motions,  To  do  this, 
it  will  be  necessary  to  make  extensive  measurements  of  hydro¬ 
dynamio  roll  moment  as  a  function  of  roll  angle,  attack  angle, 
roll  angular  velooity,  pnd  other  motion  parameters.  It  is 
recommended  that  these  measurements  be  made  using  a  rotating 
arm  facility,  which  has  been  shown  to  be  effective  for  suoh 
tests, 8®»®0»31  At  the  same  time,  the  hydrodynamic  static  and 
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damping  derivative  coefficients  should  be  measured.  Thin  will 
permit  the  values  used  for  these  coefficients  in  the  present 
study  to  be  compared  with  a  set  of  values  derived  uniformly 
from  the  results  of  a  single  static  experiment. 
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Nomenclature  for  Treating  the  Motion 
of  a  Submerged  Body  Through  a  Fluid 

Introduction 


The  itudy  of  the  motion  of  ■  rigid  body  through 
•  fluid  it  contideml  u  •  general  problem  in  the 
motion  of  rigid  bodice  under  the  action  of  gravity 
and  hydrodynamic  forces.  Thut  a  choice  of 
nomenclature  for  this  field  mutt  borrow  from  both 
rigid  dynamics  and  fluid  mechanic*, 

Like  moat  problems  in  hydrodynamics,  an  exact 
theoretical  solution  of  the  equations  of  motion  is 
available  only  for  an  ideal  fluid.  However,  as  is 
well  known,  fur  small  deviatkais  from  uniform 
linear  motion,  the  equations  of  motion  in  a  real 
fluid  can  be  linearised  and  solved  in  terms  of  a  set 
of  coefficients,  which  are  determined  either  by 
theory  or  by  experiment. 

Frequently,  rather  than  to  determine  the  trajec¬ 
tory  of  a  body,  it  suftices  to  obtain  a  measure  of 
its  stability.  A  body  in  a  state  of  equilibrium  is 


said  to  be  stable  if,  when  slightly  disturbed  from 
this  state,  it  returns  to  it,  with  perhaps  a  slight 
change  in  position  of  path.  Since  there  are  various 
possible  states  of  equilibrium,  there  are,  corre¬ 
spondingly,  various  possible  kinds  of  stability; 
these  are  defined  In  a  subsequent  section  of  this 
Bulletin. 

Although  the  motion  can  be  described  ade¬ 
quately  only  with  reference  to  some  fixed  co¬ 
ordinate  system,  it  is  shown  in  texts  on  dynamics 
that  the  equations  of  motion  of  a  rigid  body  are 
expressed  most  conveniently  in  terms  of  a  rec¬ 
tangular  coordinate  system  moving  with  the  body, 
the  so  cnllea  body  axes.  Consequently,  it  has 
been  necessary  to  introduce  nomenclature  for  de¬ 
scribing  the  orientation  of  the  body  axes  relative 
to  the  fixed  axrs. 


Body  Characteristics 


It  is  supposed  that  the  body  is  elongated  and 
has  well-defined  forward  and  after  ends,  such  as 
nose  and  tail,  or  bow  and  stern.  It  is  assumed 
also  that  the  body  has  a  “princi|xil  plane  of  sym¬ 
metry,"  as  shown  in  Fig.  I,  and  that  there  is  a 
normal  attitude  of  the  body,  with  the  principal 
plane  of  symmetry  vertical,  “o  that  the  top  (deck) 
and  bottom  of  the  body  ea  i  be  defined.  The 
body  may  be  equipped  wMi  wings,  stabilising 
and  control  surfaces,  an,’,  it  may  be  self-propelled 
or  towed. 

Srncinc  Points  of  Boot 

CB _  center  of  buoyancy  of  the  body 

CG . center  of  mass  of  the  body 

CM ,  metacenter  of  the  body 

CP . center  of  pressure  of  the  body 

CS  ,  static  center;  center  of  resultant  ol 
weight  and  buoyancy 

0.  .  ,  .origin  of  body  axes  (to  be  defined) 

TP  _ towing  point  or  towed  point 


Body  Axes 

The  origin  0  of  a  right-handed,  rectangular, 
coordinate  »y;t-.<ii  fixed  in  the  body  ia  choten 
either  to  coincide  with  CG,  when  CG  ii  in  the 
principal  plane  of  xymmetry,  or  at  any  convenient 
point  in  the  principal  plane  of  xymmetry  when  CG 
does  not  lie  in  this  plane.  Body  axes  x,  y,  s,  coinci¬ 
dent  with  the  principal  axes  of  inertia  at  0,  are 
then  defined  as  follows: 

x  .  the  kmgitudinal  axis,  directed  from  the 
after  to  the  forward  end  of  the  body 
y. .  the  transverse  axis,  directed  to  starboard 
s..  the  normal  axis,  directed  from  top  to  bot¬ 
tom  (deck  to  keel) 

The  sx- plane  is  the  principal  plane  of  symmetry. 

CoOXDINATt:»  AND  UlSTANCItS 

x„  y„  i,  coordinates  of  CB  relative  to  body  axes 
y«,  *o .  ■  coordinates  of  CG  relative  to  body  axes 
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Fio  I.— SxtrcM  of  an  Undmwats*  Boov  Bhow i No  Plann  of  Bvmmtrv  anp  Body  Ami 


x„  y»  *<  coordinater  of  CS  relative  to  body  axe* 
xTl  yf,  ar.  .coordinate*  of  TP  relative  to  body  axe* 
botB  —  the  beam  of  the  body;  a  breadth, 
width,  or  ipan 

dot  D... , diameter  of  a  body  of  revolution 


c . chord  length  of  a  win*,  rudder,  elevator, 

or  diving  plane 

m  or  GM.  .the  metacentric  height;  the  diitaoc* 
from  CG  to  CM,  positive  upward 

// . the  draft  of  a  floating  body 

k . the  depth  of  Mibmergence  of  a  tub- 

merged  body 

Jk. . . . . the  hull  height  of  a  submerged  body, 

measured  from  the  bottom  to  the  top 
of  the  hull  (pronounced  height) 

k . radiuK  of  gyration 

l  or  L . length  of  the  body 

R, . radiu*  of  turning  circle,  or  radiu*  of 

curvature  of  path 

/ . thicknew 

i  . thickne**  of  boundary  layer 


Ama*  and  Vqivmu 

A . an  area,  projected  or.  a  given  plane  or 

on  the  median  riant  of  a  hydrofoil, 
•uch  m: 

At . .projected  are*  of  bow  plane* 

A,  ■ . projected  ru  of  »tem  plane* 

A, . projected  rudder  area 

A„  A „  A, .urea  of  the  projection  of  the  tubmerged 
part  of  the  body  on  the  y»,  it,  and  xy 
plane*,  respectively 

5 . wetted  surface  area  of  the  body 

V  volume  o'  the  body  (pronounced  vol) 

Inbbtia  Characteristic* 
m . maw  of  body 

I„  I„  l,. . ,  moments  of  inertia  of  the  body  about 
x,  y,  $  axes,  respectively 

k„  k„  k,  radii  of  gyration  of  the  body  about 
x,  y,  $  axes,  respectively 


Fixed  Axes — Orientation  and  Motion  of  Body 


Fixed  Axes 

It  will  be  assumed  that  the  accelerations  of  a 
point  on  the  surface  of  the  earth  can  be  neglected 
so  that  a  set  of  axes  fixed  relative  to  the  earth  may 
be  considered  as  fixed  axes.  Choose  a  right- 
handed,  orthogonal  set  of  "fixed  axes,"  x*  y«,  *>, 
fixed  relative  to  the  earth,  so  that  the  x.  and  yo  axes 
art  in  a  horisontal  plane,  and  the  to-aris  is  vertical 
and  directed  downwards. 

x, —  the  fixed  longitudinal  axis,  in  a  fixed  direc¬ 
tion  in  a  horisontal  plane,  considered  as 
the  forward  direction 

y*. .  .the fixed  transverse  axis,  perpendicular  to 
in  a  horisontal  plane,  directed  to  starboard 
s< — the  vertical  axis,  directed  downwards 

Orientation  or  Body  Axes 

Alt  angular  displacements  will  be  taken  as  posi¬ 
tive  in  the  sense  of  rotation  of  a  right-hand  screw 
advancing  in  the  positive  direction  of  the  axis  of 
rotation  (right-hand  screw  rule). 

The  orientation  of  the  body  axes  relative  to  the 
Axed  axes  can  be  described  in  various  ways.  In 
aerodynamics  the  orientation  of  a  body  in  space  is 
usually  described  in  terms  of  an  angle  of  pitch  or 
trim  t,  an  angle  of  yaw  and  an  angle  of  roll  or 
heel  *.  In  addition  it  is  necessary  to  introduce  a 
second  set  of  angles  a,  0,  y.  which  describes  the 
orientation  of  the  body  relative  to  its  direction  of 
motion. 

Space  Orientation  8,  i,  a  and  Corresponding 
Angular  Velocities  op  Body 

The  space  orientation  of  the  body  axes  x,  y,  s 
relative  to  the  fixed  axes  x,,  y»,  to  may  be  described 
by  the  following  procedure  (see  Fig.  2).  First 
suppose  the  axes  x,  y,  s,  to  coincide  with  the  axes 
x«,  y«,  Si.  Rotate  the  body  about  s«  through  an 
angle  of  yaw  4  so  that  the  axes  Xt>,  y,  assume  the 
intermediate  positions  Xi,  yi ;  then  rotate  the  body 
about  the  new  position  of  the  y-axis  through  an 
angle  of  trim  #,  so  that  s«  moves  to  Si  and  xi  moves 
to  x;  finally  rotate  the  body  about  the  new  posi¬ 
tion  of  the  x-axis  through  an  angle  of  roll  «  so 
that  the  axes  yt,  z\  assume  their  final  positions  y, 


s.  In  accordance  with  this  procedure  we  have  the 
following  definitions: 

8. ...  the  angle  of  trim  (or  pitch) ;  the  angle  of  ele¬ 
vation  of  the  x-axis;  i.e,,  the  angle  between 
O.rand  the  horisontal  plane  xoys,  positive 
in  the  sense  of  rotation  from  the  so  to  the 
x-axis 


I,  T*  Jtsrtoard 


Hio  2.—  Oaixntation  or  B.idy  Axes  Rilativb  to 
Fixed  Axes  in  Teems  or  S,  »  #.  Viewed  peom  Below 
*y  Plane 


+... .  the  angle  of  yaw;  the  angle  from  the  vertical 
plane  Stx,  to  the  vertical  plane  to*,  positive 
in  the  positive  sense  of  rotation  about  the 
it-axis 

*. .  an  angh  of  roll;  the  angle  from  the  vertical 
xso-ptune  to  the  principal  plane  of  sym¬ 
metry  xi,  positive  in  the  positive  sense  of 
rotation  about  the  x-axis 

The  direction  cosines  of  x,  y,  s  relative  to  x,,  y«, 
to  may  be  tabulated  as  in  Table  I . 


Taui.k  1.— Direction  Cosines  of  Body  Axes  Relative  to  Fixed  Axes  in  Teems  op  0,  4,  * 


Xo 

- - X—  - 

cos  8  cos  4 

y* 

cos  8  sin  4 

So 

-sin  8 

-cos  4  sin  4  4  sin  8  sin  4  cos  4 
cos  4  cm  4  4-sin  8  sin  4  sin  4 
cos  8  sin  4 


sin  4  sin  4  4  sin  8  cos  4  cos  4 
-sin  4co»4  4*in  8  cos  4  sin  4 
cos  8  cos  4 
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Alto,  for  the  direction  contact  o f  *•  relative  to 
the  body  axet,  put  n,  -  -tin  t,  *,  »  cot  t  tin  4, 
n,  -  cot  I  cot  4, 

Let  p,  },  r  be  the  component*  of  the  angular 
velocity  vector  relative  to  the  body  Met  x,  y,  >. 
Denote  the  time  derivative*  of  t,  4,  4  by  i,  i,  4 
We  have  the  (oUowini  relatkmi: 

P  -  4  -  4  tin  # 

9  ■  4  cot  #  tin  9  +  4  cot  4 
r  -  4  cot  t  cot  4  -  4  tin  4 
p,  q,  r. .  .angular  velocity  component*  relative  to 
body  axet  x,  y,  1;  angular  velocitiei  of 
roll,  pitch,  and  yaw 

Velocity  or  Orioin  or  Body  Axat—AwoLt*  or 
Attack  a,  and  Dairr  P 

It  thould  be  recalled  that  the  origin  of  body 
axet  0  it  taken  to  coincide  with  the  center  of 
mat*  CU,  when  the  CG  lie*  in  the  principal  plane 
of  tymmetiy. 

J<et  N,  r,  w  be  the  componeuti  in  the  x,  y,  1  co¬ 
ordinate  tyitem  of  the  linear  velocity  of  0  relative 
to  the  liv'd. 

w,  r,  tv. .  component*  along  body  axe*  of  velocity  of 
origin  of  body  axet  relative  to  fluid 
U  or  V..  velocity  of  origin  of  body  axet  relative  to 
fluid,  In  general,  the  lymbol  U  it  uted 
in  this  Bulletin 

Some  of  the  hydrodynamic  force*  which  act  on 
a  body,  the  to-called  italic  force*,  depend  on  the 
orientation  of  the  body  with  respect  to  the  relative 


Pm.  3.-0*iintatiun  nr  Body  Axtt  Rhi.ativ*  to 
Velocity  or  OaiuiK 


fluid  velocity,  This  orientation  it  specified  by 
the  velocity  component*  a,  t>,  w  but  it  more  usually 
given  in  terms  of  angles  a  and  8  which  may  be 
defined  as  follows  (see  Fig,  It)  :■ 

a. . .  .the  angle  of  attack;  the  angle  to  the  longi¬ 
tudinal  body  axis  from  the  projection  into 
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the  principal  plane  of  symmetry  of  the 
velocity  of  the  origin  of  the  body  axe* 
relative  to  the  fluid,  positive  in  the  posi¬ 
tive  tense  of  rotation  about  the  y-axi* 
fi  ...  the  drift  or  sideslip  angle;  the  angle  to  the 
principal  plane  of  symmetry  from  the 
velocity  of  the  origin  of  the  body  axe* 
relative  to  the  fluid,  positive  in  the  posi¬ 
tive  tenie  of  rotation  about  the  s-axis 

We  have  the  following  relation*  between  the 
velocity  veciur  and  these  angle*  (lee  Fig.  3): 

tan  o  -  +u»/n 
•in  fl  -  —*>/£/ 

where  U  is  the  velocity  of  the  origin  of  the  body 
axes  relative  to  the  fluid. 

U*  -  N*  +  +  w’ 

Conversely,  we  have 

n  -  U  co*  a  co*  8 
r  -  -  V  sin  ft 
w  -  U  tin  or  co*  & 

Motion  Oriintation  a,  B,  y  and 
CoRRRsroNDiNO  Anoular  Vblocitik*  or  Body 

Far  bodies  towed  in  a  horixontal  straight  line, 
the  fixed  xraxit  may  be  chosen  to  coincide  with 
the  direction  of  motion.  The  orientation  of  the 
body  axet  x,  y,  1  relative  to  the  fixed  axet  xt,  yi,  * 
may  then  be  expressed  in  terms  of  angle*  a,  B,  and 
y,  defined  a*  follows: 

a  —  the  angle  of  attack,  defined  in  the  previous 
lection 

B — the  drift  or  sideslip  angle,  defined  in  the  pre¬ 
vious  lection 

T  —  an  angle  of  roll;  the  angular  displacement 
about  the  x,-*xis  of  the  principal  plane  of 
symmetry  from  the  vertical,  positive  in 
the  positive  tense  of  rotation  about  the 
xraxit 

When  the  orientation  is  specified  in  term*  of  a, 
B,  y,  the  angles  of  attack  and  drift  are  independent 
of  the  angle  of  roll.  This  it  not  the  cue  for  the 
#,  4>  4  tyitem.  Indeed  t*  -  a  and  4  -  B  only 
when  4  -  0.  For  this  reason  it  is  preferable  to 
use  the  angle*  a,  6,  y  in  cases  where  tome  i>oint 
in  the  principal  plane  of  symmetry  move*  in  a 
fixed  direction. 

The  direction  cosines  of  x,  y,  *  relative  to 
x*  yt,  u  may  be  tabulated  as  in  Table  2. 

Denote  the  time  derivative*  of  o,  8,  y  by  4,  8,  y. 
We  have  the  following  relation  between  p,  q,  r 
and  4,  f,  y  :■ 

p  -  8  sin  a  +  y  cot  a  co*  8 


Tablr  3.— Dibbction  CiwtNM  or  Boov  Axm 

- - * - >  • 

X,  cut  a  CM  0 

y*  tin  a  (in  y  4-cm  a  tin  fi  cm  y 

ti  -tin  n  cm  y  +cot  a  tin  tf  tin  y 


f  -  A  -  4  tin  fi 
r  •  4  cm  a  4-  4  tin  o  cot  fl 

Control  Surfacri 

I . angular  ditpitccmtnt  of  a  control  turftcc 

I, — rudder  ang'i,  potitlve  in  the  tente  of  poti- 
tive  rotation  about  the  t-axit 


ItRLATlVR  TO  PtXRD  AXBt  IN  TBaMB  OP  a,  6,  f 

- y - - - -* - ■ 

-tin/)  tin  a  cm  fi 

cm  f  cot  t  -cm  a  tin  y  4-  tin  •  tin  fi  cm  y 

cm  d  tin  y  cm  o  cm  t  4-  tin  n  tin  d  tin  t 


I*  angular  displacement  of  the  bow  plane,  poei- 
tive  in  the  tente  of  poeitlve  rotation 
•  about  the  y-axla 

I, —  angular  dieplacement  of  the  item  plane, 
poeitive  in  the  tente  cf  poeitive  rotation 
about  the  y-axi* 


Forces  and  Moments 


It  it  convenient  to  dittinguith  between  the 
hydrodynamic  forcet,  which  include  the  propeller 
throat  when  the  body  it  telf-propelied,  and  other 
external  forcet  auch  u  thoee  due  to  gravity 
(weight  and  buoyancy)  and  towlinet. 

At  it  the  cate  for  angular  ditpiacementt,  mo¬ 
menta  are  taken  at  poeitive  in  accordance  with 
the  right-hand  acrew  rule. 

X,  K,  Z..,. .  hydrodynanic  force  component*  rela¬ 
tive  to  body  axet,  referred  to  at 
longitudinal,  lateral,  and  normal 
forcet,  rcapcclively 

Dor  h . drag  component  of  hydrodynamic 

force  in  direction  of  relative  flow; 
may  be  qualified  by  an  appropriate 
aubKript  to  Indicate  a  component 
of  the  drag;  retiatance 

C . crot*  force;  component  of  hydro¬ 


dynamic  force  normal  to  lift  and 
drag,  poaitive  to  t tar  board 

L . lift;  component  of  hydrodynamic 

force  in  the  principal  plane  of  tym- 
metry  normal  to  relative  flow,  poei¬ 
tive  when  directed  from  bottom  to 
top  of  body  (keel  to  deck) 

K,  M,  N.  hydrodynamic  moment  component* 


relative  to  body  axet,  referred  to  at 
rolling,  pitching,  and  yawing  mo¬ 
menta,  reapectively 

Q . torque  about  the  axil  of  a  control  tur- 

face,  with  appropriate  tuOtcript 

?' . towline  teniion  (pronounced  toll) 

T„  F„  T,  .component!  of  towline  teniion  rela¬ 
tive  to  body  axe* 

W or  A. .  weight  of  the  body;  W  -  mg 
B . buoyancy  force 


W  -  B — reaultant  of  weight  and  buoyancy 


Dimensionless  Forms 


A  given  phyaical  quantity  will  be  non-dinten- 
•ion allied  by  contidering  it*  dimensions  in  term*  of 
maaa,  length,  and  time  at  fundamental  unit*,  and 
dividing  the  mate  by  (H)W,  the  length  by  I,  and 
the  time  by  l/U.  For  example,  the  dimension*  of 
ditf/dg  art  maaa  X  (length)1  X  (time)-'..  Hence 

dA/ _ t _ l  f  dM/tXj 

*  (HWf*  V  ”  mwv 


I  m  it 


*• 

1 


When  a  itability  analytia  hat  been  completely  ex- 
pretted  in  non-dimentional  form,  the  prime*  may 
be  omitted,  provided  that  a  itatement  to  that 
effect  it  made. 


Bodv  Shapr  Parambtrrx 


it  the  corresponding  dimenaionlttt  exprettion. 

The  non-dimentionalUed  form  of  a  given  phyti- 
cal  quantity  will  be  indicated  by  the  prime  of  that 
quantity,  unlett  explicitly  defined  otherwiae. 
Thu* 


The  thape  of  a  body  may  be  deacribed  conveni¬ 
ently  in  termi  of  the  ratio#  of  tome  of  its  principal 
dimension*,  such  at: 


A -7 


b'  -  b/l 


beam-length  ratio 


<Y m  d/l . diameter-length  ratio  of  a  body 

of  revolution 

1/4 . flneneaa  ratio 

IV  «  H/I . draft-length  ratio 

b/H . beam-draft  ratio 


C,  ••  A, /bit  .  .maximum  Motion  coafHdant 
S'  -  S/I'. ..... .a  wctted-surface  coefficient 

S/Wi.  a  wetted-surfacc  coefficient 
V  -  V/l'. , .  volumetric  coeflicient 
C,  -  V HA,  ,  priimatic  coefficient 
C,  -  V/MI — block  coefficient 


u  -  b'/A, . aspect  ratio  of  a  wing 

Vrlocitira  and  Accrlrrationr 

a'  -  u/U,v'  -  v/U, 

w'  -  w/U . dimensionless  velocity 

components 

*'  -  *1/1/',  t'  -  H/U', 
w'  •  sif/f/1 .  dimensionless  accelera¬ 

tion  components 

P'  -  pf/(/.  «'  -  gl/f/, 

r'  •  rl/U  . dimensionless  angular 

velocity  components 

- 

I'mH'iV*  .  .  ..dimensionless  angular 

acceleration  com¬ 
ponents 


Forcr  and  Momrnt  Corfmciknta 

Since  forces  and  moments  are  frequently  non- 
dimenaionalised  in  terms  of  S,  A,  or  I',  alternate 
forms  will  be  presented. 

X'  »  X/(\  ilfl'U'  kmgitudinul  force  coefficient 

V  m  Y/(l/i)pl'U'.  lateral  force  coefficient 
Z‘  m  Z/('/i)pl'V '  .  normal  force  coefficient 

D‘  •  D/('/Z)pl'U'.. ,  .drug  or  resistance  coefficient 
Ct,  «  D/('/i)pAU'  .  .drag  or  resistance  coeffi¬ 
cient,  alternate  form 

to  "  D/(\%)pY','U'  .Ata%  of  resistance  coeffi¬ 
cient,  alternate  form 

C,  -  D/0/£)pSU' — drag  or  resistance  coeffi¬ 
cient,  alternate  form 

C /  -  D//('/i)i>SU'  friction  drag  or  resistance 
coefficient 

C,  -  D,/('/i)pSU'  .residuary  drag  coefficient; 

C  ~  C  —  C 

C'  »  C/dftpl'U', . .  cross  force  coefficient 
Cr  m  C/M)pA  U'.  cross  force  coefficient,  alter¬ 
nate  form 

V  -  L/m)pl'U'.  lift  coefficient 

Ci  -  L/\'.i)pAU'.  lift  coefficient,  alternate 

form 

K’  -  Ki(\:i)pI'V' .  . rolling  moment  coefficient 
C*  “  A.'/C  ■i)pAlV'  rolling  moment  coefficient. 

alternate  form 

.1/’ «  A//(*  i)pl'U'  pitching  moment  coefficient 


Cj»  -  M /('  typAW'.  pitching  moment  coeffi¬ 
cient,  alternate  form 

N'  -  N/{A)M*U'.  .  yawing  moment  coefficient 
C,  -  N/(K\\tAW'. .  .yawing  moment  coefficient, 
alternate  form 

CORFFtCIRNTR  OR  STATIC  FDRCR  DRniVATIVR* 

The  partial  derivative  of  a  force  or  moment  com¬ 
ponent  X,  Y,  Z,  K,  M,  or  N  with  respect  to  a 
linear  or  angular  velocity  or  acceleration  a,  v,  w, 
P,  ?,  r,  rt,  #,  W,  fi,  f ,  t  will  be  designated  by  the  force 
or  moment  with  the  velocity  or  acceleration  as  a 
subscript;  e.g.,  HY/bp  -  V,. 

Y  1  m  —*  _  V'  _  y  >  „  Xw 

*  (JflgW'*'  ('.$&>'• 

V  I  m  1 »  _  J’ '  m  „  Yi Y  I  m  jj» 

*  OgWI/*  OAW'LI'  "  (Wi/ 


•  i 

,, 


2 

**  “  (Rw*i/‘ z:  w  (bj^*  •'  "  (%yw 

Note  that  also  X.'  -  dX'/du\  X.'  -  dX'/fr>\ 
etc. 


CoRrriciRNTS  or  Kotarv  Porcr  Drrivativrs 

-Is 

~  c w'v 


Y  1  _  'It  -  Y  1  _  -Jis-  V  '  _ 


IHWl' 

t-  >  «,  I’l  „  !i  V  '  m  ~  2j__ 

7  '  _ _ “7  .  X,  _ ,  z, 

'  OiW  '  (by't"  (tg’yt/ 

Note  that  also  .V,'  -  dA'/dp',  X,'  -  &A"/<V, 
rtc. 

V IRTUAI.  INKRTIA  CoKPriCIRNTS 


-V.' 

-  ,  X> 

CiV* 

-  -Y'  Y  '( 

0*)^  * 

X* 

(b)W* 

»v 

_  b  y ' 
(bW  * 

-  *  '  ,  K*' 

_  r* 

(bW 

z: 

2*  ,, 
“  Ciy*  ' 

_  »  f 

;bw*’ 

- 

053* 

xt‘ 

"  (WW*1  A* 

-  Y'  V,' 
CvV 

x, 

(J<W* 

»v 

„  1*  ... 

“  , ,JL-.  IV  - 

J’r 

(b)>*  * 

CiW4 

*  (4W* 

„  2>  „ , 
C  Ypf  * 

-  -?*  z, 
CiW'  ' 

"  (b)> 

Note  thut  also  AY  » 

d.V'  dii’,.  XV 

-  dX'  dp', 

etc. 


Corfficirnt*  of  Static  Momrnt  DrrivAtivks 


Note  that  alio  K,'  -  bK'/du',  K,  -  HK'/d 
etc. 

Virtual  Momrnt  or  Inrrtia  Corfficirntr 


Note  that  alio  K -  DK'/W,  K,'  -  dK'/dp', 
etc. 

CORFFICIBNTR  OF  ROTARY  MoMRNT  DrRIVATIVRR 


«  .  ,  O'* 

(HW4’  *  (H)af4 


Note  that  alio  K,'  -  d/C7d/>',  X,' 
etc. 


d*7*\ 


Linearized  Expansion  of  Force  and  Moment  Coefficients 


In  the  following  analysis,  all  quantities  are  non- 
dimeniionallied,  «o  that  primes  may  be  omitted. 

It  will  be  assumed  that  the  hydrodynamic  force 
and  moment  coefficients  depend  only  upon  the 
non-dimensionalieed  linear  and  angular  velocity 
and  acceleration  components;  i.e.,  X  -  X  (u,  t>, 
*>,  p,  9>  r>  *■  «*>  P>  ii  0.  with  similar  forms  for 
Y,  Z,  K,  M,  and  N. 

Let  u„  *„  w„  p„  q„  r„  ti„  *„  re,,  p„  $„  f,  describe 
the  state  of  motion  of  the  body  at  s  time 
Then,  for  small  changes  from  this  initial  state  we 


have  the  linearised  Taylor  expansion 

X  ■  X,  +  Xm(u  -  u.)  +  *„(»  -  y.)  + 

*„(«>  -  re.) 

+  *„(/>  -  P.)  +  XJq  -  ,.)  + 

AT„(r  -  r.) 

-I-  -  d.)  +  X„(t  -  *,)  + 

-  *,) 

+  **(/>  -  />*)  +  X*.(i  -  4t)  + 

*  X,.(>  -  K) 

with  similar  expressions  for  Y,  Z,  K,  M,  and  N. 


Equations  of  Motion  for  a  Free  Body 


The  simplest  form  of  the  equations  of  motion  is 
obtained  with  body  axes  coincident  with  the  prin¬ 
cipal  axes  of  inertia,  and  the  origin  ut  the  renter 
of  mass  CC.  For  this  case  the  equations  are,  in 
dimensionless  form, 

X  -  !»[<J  +  qw  -  rv| 

Y  -  mj*  +  ru  —  pw\ 

Z  -  ih|w  +  pv  -  ?u| 

K  -  l,p  4-  (/,  -  l,)qr 
Af  m  ij  +  (l,  -  l,)rp 
N  “  l,r  +  (/,  --  f,)pq 

When  the  origin  of  the  body  axes  is  not  at  CO, 


as  has  been  |>roposed  when  CG  does  not  lie  in  the 
plane  of  symmetry,  the  equations  of  motion  as¬ 
sume  the  more  complex  form  . 

.V  -  m|d  •+■  qw  -  rv  -  xa(q •  -4-  r*)  + 

yjpq  -  f)  +  *dpr  +  i)l 
with  similar  expressions  for  Y  and  Z. 

X  -  l,p  -i-  (/,  -  l,)qr  + 

+  pv  -  qu)  -  *,,($  +  ru  -  pw) | 

with  similar  expressions  for  M  and  N. 

The  linearization  of  these  equations  in  a  particu¬ 
lar  problem  introduces  considerable  simplification. 
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Equilibrium  Equations  for  a  Towed  Body 


The  application  at  the  propoeed  nomenclature 
to  the  caae  at  a  (table  kite  or  paravane  will  be 
considered  aa  an  example.  Equation!  (or  deter¬ 
mining  the  attitude  ot  the  body  when  the  hydro* 
dynamic  characteristics  ot  the  body  arc  known 
will  be  presented. 

It  will  be  supposed  that  the  body  is  towed  at  a 
uniform  speed  in  a  Used  direction  In  the  horisontal 
plane.  Let  the  xraxis  be  In  the  direction  ot  tow. 

Since  equilibrium  is  assumed,  the  hydrodynamic 
forces  and  moments  are  functions  ot  m,  t>,  and  w, 
or  a,  4,  U  only;  e.g.,  X  -  X  (»,  p,  re). 

The  resultant  of  weight  end  buoyancy  is  the 
force  W  -  B  in  the  at  direction,  acting  at  C5, 
The  direction  cosines  of  the  Sraxis  relative  to  the 
body  axis,  n„  n„  *„  in  terms  of  a,  4,  y,  may  be 
rtad  from  Table  2. 

The  towlinc  tension  sets  ut  a  point  TP  whose 
coordinates  (*,,  >r,  s f)  are  functions  of  a,  4,  y 
which  depend  upon  the  type  ot  bridle  used  for 
attaching  the  towline  to  the  body,.  In  the  sim¬ 


plest  case,  when  the  towline  is  attached  directly 
to  the  body  (single  point  of  attachment),  we  may 
take  the  origin  ot  body  axes,  0,  at  TP,  so  that 

*r  -  yf  -  *r  -  0. 

The  equations  may  be  written  as  follows: 

*(«,  4,  V)+T,+  n,(W  -  B)  -  0 
V(«,  4,  U)  +  F,  +  n,(lt'  -  B)  -  0 
B(< »,  4,  U)+T,  +  ».(IV  -  B)  -  0 

K(a,  4,  U)  +  y,T,  -  srF,  + 

(IV  -  B)(y,n,  -  s,*»,)  -  0 
!*(«,  4,  U)  +  s,F,  -  xrT,  + 

(IV  -  B)(s.x.  -  x,n,)  -  0 
N(«,  4,  V)  xfT,  -  yrF,  + 

(IV  -  B)(x,n,  -  y,H.)  -  0 

Since  n„  m„  n,  are  functions  of  a,  4,  y  from  Table 
2,  and  xTl  y„  sr  are  also  assumed  to  be  known 
functions  ot  a,  4,  V  from  the  bridle  geometry, 
these  are  six  equations  In  the  unknowns  a,  4,  y> 
F„  F„  F,. 


Physical  and  Other  Terms 


a.- . mass  density 

g .  acceleration  due  to  gravity 

F  m  U/Vtl,  V/Vgi.Froude  number 

a . coefficient  of  viscosity 

r . kinematic  viscosity,  p/p 


X) 

X. 

R. 


vi/*,  vi/». 

Ud/r,  Vd/r. 
Ux/>,  Vx/¥., 
f/l/r,  VI/*, 


Keynolds  numbers  corre¬ 
sponding  to  various  linear 
dimensions 


X . linear  ratio  of  full  site  body 

to  model 

/ . time 

a . a  frequency 

S  -  hI/U .  Strouhal  numbers  c«>rre- 

S,  -  tui/U  sponding  to  various  linear 

dimensions 

» . cavitation  number 

. . roots  of  the  stability  equa¬ 
tion,  *  »  I,  2,  .  .  . 


Definitions  of  Stability 


Mktackntbic  Stahiuty 

A  body,  floating  in  equilibrium  either  com¬ 
pletely  submerged  or  on  the  surface  of  a  fluid,  is 
said  to  have  melactHlrU  stability  if,  when  disturbed 
from  the  equilibrium  position  In  either  trim  or 
heel,  it  returns  to  the  original  pwition.  If  it  is 
stable  aguinst  a  disturbance  in  heel,  it  is  said  to 
have  tranmrse  mttactntrir  stability.  If  it  is  stable 
against  a  disturbance  in  trim,  it  is  said  to  have 
lonnilutlinal  mrlairntrU  stability.  A  body  may 
have  metneentrio  stability  about  a  position  other 
than  the  upright  position;  in  thi;.  case,  the  posi¬ 
tion  of  equilibrium  must  lx1  stuted.  Metacentric 


stability  implies  that  the  metacentric  height  m 
is  positive;  that  is,  CM  is  above  CG. 

Tbanslational  Stability 

In  discussing  the  directional  stability  of  a  body 
it  is  convenient  to  distinguish  between  the  stability 
of  the  body  when  restrained  to  move  only  like  a 
weathercock,  and  the  stability  of  the  unrestrained 

motion. 

(a)  StatU  ( or  Weathercock)  Stability.  Suppose 
the  body  to  be  moving  in  a  constant  direction  and 
restrained  so  that  its  only  freedom  of  motion  is 
that  of  rotation  about  an  axis  perpendicular  to 


A-iO 


the  direction  of  motion.  The  body  to  Mid  to  be 
•Utterly  itebto  about  the  given  axle  in  eome 
equilibrium  orientation  about  it,  if  it  return!  to 
thle  orientation  after  a  (light  disturbance.  If  the 
body  to  (table  about  the  transverse  axto,  it  to  Mid 
to  have  static  pitching  stability;  in  thto  caM,  the 
dope  of  the  pitching  moment  curve  M,  to  negative. 
If  (table  about  the  normal  axto,  it  to  Mid  to  have 
italic  yawing  stability ;  In  thto  caxe,  the  (lone  of 
the  yawing  moment  curve  Nt  to  negative.  Static 
(tability  require!  that  the  center  of  preeeure  be 
aft  of  the  axto  of  rotation.  The  negative  magni¬ 
tude!  of  the  derivative!  Mt  and  N,  for  aome  oriett- 
tation  of  a  body  are  a  meaetire  of  the  static  (tabil¬ 
ity  of  the  body.. 


(*)  Dynamic  Stability  an  Cantu.  A  body  to 
Mid  to  be  dynamically  stable  an  cantu  if,  after  its 
steady  state  motion  hat  been  (lightly  disturbed 
from  a  straight  count  with  Axed  control  surfaces, 
it  resume!  iu  motion  on  another  straight  course. 

\  (hip  to  said  to  be  dynamically  stable  in  a  turn 
If,  after  its  motion  hat  been  (lightly  disturbed 
front  a  given  circular  turn  with  Axed  control 
surfaces,  it  resumes  a  motion  in  a  turn  of  the 
Mme  diameter,  with  a  slight  displacement  of  the 
path. 

Dynamic  stability  can  be  further  claseifled  aa 
longitudinal  stability,  for  motion  in  the  plane  of 
symmetry;  or  kHatal  stability,  lot  motions  involv¬ 
ing  rolling,  yawing,  and  sideslipping. 


Summary 


Operations 
(dot)  i.ver  t  symbol 
'  (prime)  of  a  symbol 
subscript  u, »,  w,  p,  f,  r,  X,  S,  ti>,  p,  i  or 
I  applied  to  X,  Y,  Z,  K,  At,  nr  N 
'  (prime)  of  result  nf  previous  opera¬ 
tion 


Typical  non-dimensional 

Symbol 

formula 

a  -  b'/A, 

A 

A'  -  A//* 

At,  A„  At 

At  -  At/I' 

A#,  Af,  A, 

A,'  -  A, /I' 

C,  -  A, /bit 

b 

y  -  bn 

B 

B'  -  B/M)fl'U' 

c 

('  -  c/l 

c 

C‘  -  C/M)a)'U> 

Cc  -  C/M)aAU' 

CB 
CG 
CM 
CP 

cs 

d  or  D 
Dor  R 


0/ 

n, 

F 


h1  -  d/i  -  n/i 

D‘  -  DKMuPV* 
CK- D/M), A  U' 
ft  -  D/M),  V'>U' 
C,  -  D/M)l>SV' 

C,  -  D,/M  ),SV' 
c,  -  n,/M),FV' 


Example  Definitions 

X  -  dn/dt  Derivative  with  respect  to  time 
«*'  ■  m/M)^1  Non-dimensional  form  of  a  symbol 

Y,  •  BY/Bp  Partial  derivative  with  respect  to  sub¬ 

script 

Y,'  -  (BY /Bp)'  Non-dimensional  form  of  the  partial 

-  BY' /Bp'  derivative  with  respect  to  the  sub¬ 
script 

Definition  of  symbol 

Aspect  ratio 

A  projected  area,  qualified  by  an  appropriate 
subanipt 

Projected  area  of  bow  plant,  rudder  and  stern 
plane,  respectively 

Area  of  projection  of  submerged  part  of  body  on 
yt,  ix,  and  xy-planes,  respectively 
Hearn  of  body 
Buoyancy  force 
Chord  length 
Cross  force 

Center  of  buoyancy 
Center  of  mass  of  body 
Metacenter  of  body 
Center  of  pressure 

Static  center;  center  of  resultant  of  weight  and 
buoyancy 

Diameter  of  a  body  nf  revolution 
Drag  forre;  resistance 


Friction  drag 

Hesiduary  drag,  C,  -  C,  -  C, 
Proude  number,  U/Vgl,  V/\/gt 
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Symbol 

Typical  non-dimensional 
formula 

1 

. 

mot  GM 

H 

//'  -  //// 

h 

It 

l„  /, 

-  /./(H)gf* 

h 

*'  -  *// 

« *,// 

K,  M,N 

A'  -  A/Wi'* 

A'. 

C„  -  K/iHlxAlU'  , 

A,'  -  K,/{WU  m  dK‘/du 

As 

As'  -  A*/(H)pf‘  -  dK'/dtt1 

A’s 

A,'  -  A,/()  sip/4!/  -  dK'/dp' 

As 

As'  -  A, /(i  4)sf*  -  dK'/dp' 

/  or  L 

T  -  1 

l 

/-'  -  L/(h)M'V% 

m  or  GM 

Ci  -  l./(h)pAV' 

m 

m'  -  w/Oj)#/* 

n 

»'  -  nl/ir 

n„  n, 

0 

Pi  *>  ' 

p'  -  pl/U 

Q 

Q'  -  4>/(l»)sf(/' 

K  or  D 

K, 

s 

5’  -  ,V//' 

s»st 

1 

r  -  t/l 

1 

1'  -  <(//< 

r 

.r  -  r/(H)W' 

T„  r„  r, 

. 

TP 

u,  s,  w 

u'  -  «/(/ 

(/or  V 

(/'  -  1 

V 

V'  -  v//> 

Word 

C.  -  V/lbll 

Cr  -  H/M, 

W'  -  W/(\<j)pl'V' 

x,  y,  i 

*'  -  *// 

Definition  of  symbol 
Acceleration  of  gravity 
Metacentrtc  height 
Draft  of  a  body 

The  depth  of  tubmergence  of  a  submerged  body 
The  hull  height  of  a  submerged  body,  measured 
from  the  bottom  to  the  top  of  the  hull  (pro* 
flounced  height) 

foments  of  inertia  about  x,  y,  s  axes,  respectively 
Radius  of  gyration 

Radii  of  gyration  about  x,  y,  t  axes,  respectively 
Rolling,  pitching,  and  yawing  moments,  respec¬ 
tively 

Typical  static  moment  derivative;  derivative  of 
a  moment  component  with  respect  to  a  velocity 
component,  0  K/du 

Typical  moment  of  inertia  coefficient;  derivative 
of  a  moment  component  with  respect  to  an 
acceleration  component,  dK/dd  . 

Typical  rotary  moment  derivative  of  a  moment 
component  with  respect  to  an  angular  velocity 
component,  UK /dp 

.  Typical  moment  of  inertia  coefficient;  derivative 
of  a  moment  component  with  respect  to  an 
angular  acceleration  component,  dK/dp 
A  characteristic  length  of  body 
Lift  force 

Metacentric  height 
Mass  of  body 
Frequency 

Direction  cosines  of  vertical  relative  to  body  axes 
Origin  of  body  axes 

Angular  velixrltlcs  of  roll,  pitch,  and  yaw,  respec¬ 
tively 

Torque  almut  axis  of  a  control  surface 
Resistance  or  drag  force 
Radius  of  turning  path,  or  radius  of  curvature 
Reynolds  number,  t///»,  f/rf/r,  Ux/t,  Vt/» 

Wetted  surface  area 
Strouhal  number,  nl/U,  hH/U 
Thickness 
Time 

Towline  tension  (pronounced  toll) 

Components  of  towline  tension  relative  to  body 
axes 

Tow  point  at  body 

Longitudinal,  transverse,  and  normal  compo¬ 
nents,  respectively,  of  the  velocity  of  the  origin 
of  body  axes  relative  to  the  fluid 
Velocity  of  origin  of  body  axes  relative  to  fluid 
Volume  of  body  (pronounced  vol) 

Block  coefficient 
Prismatic  coefficient 
Weight  of  body 

Body  axes,  or  c<x»rdinatr»  of  a  point  relative  to 
body  axes 
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,  Symbol 

Typical  non-dimensional 
formula 

Definition  of  aymbol 

i  *»  y»  h 

. 

Coordinate*  of  center  of  buoyancy  relative  to 

«•.?*** 

*#'  -  *«// 

body  axei 

Coordinate*  of  center  of  mat*  relative  to  body  axes 

1 

Pixed  axes,  or  coordinate*  of  a  point  relativ*  to 

1 

*•.  y*  •• 

fixed  axe* 

Coordinate*  of  static  center  relative  to  body  axes 

*  *r ,  Vr,  b 

Coordinate*  of  towpoint  (at  body)  relative  to 

1  X,  Y,Z 

X'  -  */<H)*W 

body  axes 

Lonfitudinal,  lateral,  and  normal  components, 

1  *• 

*.'-*,/W(/-dA(7dn' 

respectively,  of  hydrodynamic  force  on  body 
Typical  static  fore*  derivative;  derivative  of  a 

1 

Xt 

Xi'  -  Jf./(H)W  -  Wf'/as' 

force  component  with  respect  to  a  velocity 
component,  bX/b# 

Typical  inertia  coefficient;  derivative  of  a  force 

1  , 

x,'-x,/(W[/-dX'/dp' 

component  with  respect  to  an  acceleration  com¬ 
ponent,  bX/bit 

Typical  rotary  force  derivative;  derivative  of  a 

1 

x,'  -  -  &x'W 

force  component  with  respect  to  an  angular 
velocity  component,  bX/bp 

Typical  inertia  coefficient;  derivative  of  a  force 

I  m,0, y 

component  with  respect  to  an  angular  accelera¬ 
tion  component,  bX/bp 

Angles  of  attack,  drift,  and  roll,  respectively 

■  « 

Thickness  of  boundary  layer 

< 

Angular  displacement  of  a  control  surface 

1  I,,  It,  1, 

Rudder  angle,  bow  plane  angle,  stern  plane  angle 

1  A 

Weight  of  body 

#,  #,  * 

Angles  of  pitch,  yaw,  and  roll,  respectively 

X 

Linear  ratio,  full-scale  sire  to  model  sixe 

1 

Coefficient  of  viscosity 

1 

Kinematic  viscosity,  a/a 

a 

Mass  density 

m  * 

•  Mill, 

Cavitation  number 

1 

...... 

Roots  of  stability  equation,  i  ■  1,3, 

I 


OTHER  SYMBOLS 

Other  symbols  used  in  this  report  not  inoluded  in 
the  SNAME  Technical  and  Research  Bulletin  1-5  (pp.  A-l  through 
A- 13)  follow: 

Dimensionless  Forms.  The  characteristic  area  used 
in  forming  dimensionless  quantities  is  A x,  the  maximum  oross- 
sectional  area  of  the  bare  hull.  The  length  l  is  the  overall 
length  of  the  missile. 

Origin  of  Body  Axes.  The  origin  of  the  body  axes  is 
considered  to  be  tho  center  of  buoyancy,  CB,  of  the  hull. 


A,B,C,D,A*,B*JC*  constants  appearing  in  approximate  form 

of  high  speed  vertical  planar  motion 
equations 

(A  projected  area  of  (2)  fins 

p  v 

AR1  aspect  ratio  of  fins  -  augmented  to 

account  for  fin-body  interference 
effects  when  used  in  Prandtl  formula 
for  zj. 

b  constant  appearing  in  solutions  of 

motion  equations 

bfc  tip-to-tip  fin  Bpan 


C  body  interference  vertical  force  coef- 

ZB1  ficient 


fin  normal  force  coefficient 


0  fin  normal  force  derivative  coefficient 

za  based  on  (A  )„ 

P  * 

D  .  R  constants  appearing  in  representation 

p  **  for  X£  based  on  the  drag  characteristics 

of  a  flat  plate 

D  ,  R  constants  appearing  in  representation 

£  for  X^  based  on  the  drag  characteristics 

of  a  sphere 
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moment  of  inertia  of  displaced  fluid 


constants  appearing  in  approximate  solution 
of  motion  equations  for  high  speed  vertioal 
planar  trajectories 

“v 

added  mass  of  parallelepiped  -  used  to 
estimate  virtual  masses  and  moments  of 
inertia  of  fins 

“l 

apparent  mass 

V  qp 

homogeneous  and  particular  components, 
respectively,  of  solution  for  q  obtained 
from  approximate  form  of  high  speed  vertioal 
planar  motion  equations 

•i»e«it1 

dimensions  of  parallelepiped,  in  om  -  used 
in  calculation  of  added  mass 

*o 

initial  value  of  time  in  trajectory  pre¬ 
dictions 

X' 

°b 

base  dreg  coefficient 

x» 

°f 

friction  drag  oosffioient 

x»p 

drag  coefficient  function  based  on  drag 
characteristics  of  a  flat  plate 

drag  oosffioient  function  based  on  drag 
characteristics  of  a  sphere 

( Z^)  bh» 
(Nq)jgl»  etc. 

hydrodynamic  coefficients  of  bare  hull 

^Xo^P*  ^Z^F* 

(MI) etc. 
q  P 

hydrodynamic  coefficients  of  fins 

xb 

x- coordinate  of  base  of  hull 

XP 

x- coordinate  of  point  of  application  of 
(zi)_ 

w  F 
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m 


I 


i 


i 

t 

I 

\ 

t 

1 

r 

& 

| 

s 

| 

l 


xx  x-coordinate  of  point  of  application  of  zx 

zx  increase  in  normal  foroe  derivative  due  to 

the  addition  of  fins  to  the  bare  hull 


Amxi  AX£, 

Au(o),  etc.  errors  in  values  of  input  parameters 

At  iteration  period  in  numerical  solution  of 

motion  equations 

Axq,  Azq,  etc.  errors  in  solutions  of  motion  equations 


C*  T) 


's' 


X#  0 


\,HA**»* 


0 1  °2 


T 


constants  appearing  in  straight  vertical  tra¬ 
jectory  motion  equation  with  constant  drag 
coefficient 

constants  appearing  in  straight  vertical  tra¬ 
jectory  motion  equation  with  drag  coefficient 
function  based  on  drag  characteristics  of  a 
sphere 

constants  appearing  in  solution  of  straight 
vertical  trajectory  motion  equation  with  drag 
coefficient  function  based  on  drag  character¬ 
istics  of  a  sphere 

constants  appearing  in  approximate  solution 
of  motion  equations  for  high  speed  vertical 
planar  trajectories 

roote  of  indicial  equation  in  approximate 
solution  of  motion  equations  for  high  speed 
vertical  planar  trajectories 

function  appearing  in  approximate  solution  of 
motion  equations  for  high  speed  vertical 
planar  trajectories 


R-898 

A-16 


APPENDIX  B 

EVALUATION  OP  HYDRODYNAMIC  DERIVATIVE  COEFFICIENTS 


1.  Evaluation  From  Experimental  Data 


a.  Pairing  of  Curves 

Certain  datr.  were  presented  in  the  form  of  graphe 
of  hydrodynamic  force  and  moment  coefficients  plotted  versus 
angle  of  attack.  To  evaluate  derivative  coefficients  from 
these  data,  curves  were  faired  and  the  values  of  their  slopes 
or  intercepts  were  estimated.  Since  the  curves  were  faired 
by  ’’eye, "  the  prooess  was  repeated  several  times,  by  two  dif¬ 
ferent  people.  The  reported  values  for  the  derivatives  are 
averages  of  the  repeated  evaluations. 

b.  Conversion  of  Coefficients  from 
Aerodynamic  Reference  Frame  to  x,  y,  z  System 

Forces  and  moments  measured  with  respect  to  a 
standard  aerodynamic  reference  frame  (see  sketch)  may  bo  con¬ 
verted  into  hydrodynamic  coefficients  with  respect  to  the 

x,  y,  z  frame  with  coin¬ 
cident  origin  by  applying 
the  following  three 
equations: 

X*  -  -CD  cos  ort-C^  sina, 

-  -Cp  sin  a-C^  cosa, 

M  -  -CM 

In  the  present  case,  however,  it  was  more  convenient  to  trans¬ 
form  the  derivative  coefficients  themselves.  Using  the  pre¬ 
ceding  equations,  taking  the  required  derivatives  and  evalua¬ 
ting  them  at  a  =  C,  gives  the  necessary  relationships: 


M'  =  -C 
a 


M 


a 
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0. 


Conversion  of  Derivatives  with  respect  to  a 
into  Dsrlvstlves  with  reepeot  to  w» 

From  Fig.  3  of  Appendix  A# 
w1  -  jj  ■  sin  a 
Therefore*  by.  the  oheln  rule  for  differentiation* 

A'  l*  &m  <oot  a> 

Thus*  for  a  »  0* 

Zw"  Za 

Nr-Mi 

eto.  t 

d.  Conversion  of  yaw  Plane  Derivative 
into  Pitoh  Plane  Derivative 

Because  the  xy  plane  ana  the  xe  plane  are  similar 
planes  of  symmetry  for  the  Basic  Pinner*  the  yaw  plane  and 
pitch  plane  hydrodynamic  derivatives  are  directly  related  by 
the  following  formulas*  easily  derived  by  geometric  considera¬ 
tions. 


z*  -  y» 

W  V 

N 

■ 

►< 

M'  -  N* 

> 

3r 

q  r 

z »  «  -y » 

Mi  -  Ni 

•  q  r 

xwq  "  "*vr 

q  r 

e.  Conversion  of  Derivatives  to 
_ x*  y,  z  System  with  Origin  at  OB 

Let  quantities  measured  wi/;h  respect  to  the  x*y,z 
frame  whose  origin  is  at  the  CB  be  identified  by  the  subscript 
"r."  Consider  another  x,y,z  coordinate;  system  whose  origin 
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has  coordinates  (xts,0#0)  in  the  V  frame.  Let  quantities 
measured  with  respect  to  this  frame  be  identified  by  the  sub¬ 
script  "t."  The  hydrodynamic  derivative  coefficients 
measured  with  respeot  to  the  "t"  frame  are  converted  to  the 
V  frame  by  application  of  the  following  equations.81 

<z;>.  -  <z;>t 

<zi>.  -  <zj>t 
<zi>.  ■  <zj)t  -  *i.  <z;)t 
*  K>t  -  *{. 

Ci>.  •  •  *{.  zi> t  +  +  *{.*  <zi>t 

<H1>.  ■  <"i> t  -  *{.  +  +  *{.*  <zi>t 

L 


f.  Conversion  of  Derivative  Coefficients  to 


Stand?t>d 


The  oharaoteristlo  area  used  herein  in  forming 
dlmenslonlee i  quantities  is  Ax#.  the  maximum  orose-seotional 
area  of  the  bare  hull.  The  length  K  is  the  overall  length  of 
the  missile.  Derivative  coefficients  based  on  different  non- 
dimenslonalising  parameters  are  converted  to  this  basis#  by 


multiplying  them  by  appropriate  ratios.  For  example#  if 

(MJJ).  a  denotes  the  ooeffloients  of  damping  momet.t  derivative 
q  Ax,  * 

based  on  maximum  cross-sectional  area  and  maximum  diameter# 
l.e.#  Mq  -  ^  pAxda  (M^)a  d  the  conversion  is  made  by  applying 


<  -  (f)’  in-), 


where  denotes  the  coefficient  based  on  Ax  and  K, 
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2.  Sample  Calculation 

Several  of  the  procedurea  described  In  the  preced¬ 
ing  aeotlon  now  will  be  illustrated  in  determining  one  of  the 
derivative  coefficient  values.  The  example  shows  the  deriva¬ 
tion  of  the  value  of  based  on  data  for  a  4MD  wood  model 
reported  by  Heald  and  Adams.8 

The  relevant  test  data  are  presented  in  a  graph  of 

[(M*).  vs  a,  where  the  interior  subscripts  indicate  the 

Vdlt 

nondimensionallzing  parameters,  Ax  and  d,  and  the  exterior 
subscript  t  indicates  that  the  origin  of  coordinates  of  the 
x,y,a  system  is  not  at  the  CB,  but  at  the  base  of  the  missile. * 


Repeated  fairings  of  these  data  by  two  different 
people  yielded  estimates  of  elopes  in  the  range,  23.0  <; 
(M^  ,<ilt  1  2^.l  (a  ln  radians).  The  value  used  in  the 


*  This  la  not  the  notation  of  ref.  5.  It  has  been  Introduced 
hero  for  consistency  with  the  preceding  sections. 
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subsequent  calculations  Is  |^M£)a  ^Jt  -  24.0.  This  value 
corresponds  to  the  line  shown  on  the  graph  above. 

To  convert  to  the  nondimeneional  form  based  on  Ax 
and  i,  the  rules  on  p.  B-3  are  used  to  give 

4  •  <!>  [(Miv]e  •  <*!>  <24'°)  •2-1'0 

Then,  using  the  results  on  p.  B-2, 

[«i]t  ■  [Mi]t  •  2-40 

To  adjust  this  value  to  the  x,y,s  coordinate  system 
with  origin  at  the  CB  by  the  methods  on  p.  B-2,  it  is  neces¬ 
sary  to  know  the  corresponding  value  for  the  coefficient  of 
statlo  foroe  derivative,  Z^.  Repeated  fairings  of  the  Z' 
vs  a  curve  and  application  of  the  results  on  p.  B-2  yielded 
the  range  -12.0  £  Z£  £  -11.8.  Using  the  value  Z^  -  -11.9 
together  with  x£s  -  -  0.40,  the  dimensionless  longitudinal 
coordinate  of  the  base  relative  to  the  CB  yields  the  final 
result 

[Hi],  -  ["ijt  -  *{.  [z;]t  ■  t-11*  -  (••w-ii.j)  -  -2.* 


3.  Theoretical  Estimates 


The  hydrodynamio  static  and  damping  derivative  co¬ 
efficients  were  estimated  by  means  of  relationships  based  on 
geometric  arguments.  These  equations,  similar  to  many  geo¬ 
metrically  founded  relationships  commonly  used  in  theoretical 
aerodynamics,1*  were  formulated  in  the  form  below  by  Strumpf.1 


zq  "  1 z<pBH  "  *1*1 

K  "  (Mw}BH  ■  *1*1 
+  *18*1 


( B-l) 
(B-2) 

(B-3) 

(B-4) 
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where  tha  ooaffioianti  with  aubaoripta  MBH1'  repreaent  bare 
hull  valuaa,  s[  ia  tha  inoreaae  in  tha  statio  foroa  deriva- 
tiva  ooaffioiant  dua  to  tha  addition  of  fin*  to  tha  bara 
hull i  and  x|  ia  tha  dinanaionlaaa  longitudinal  ooordlnata  of 
tha  point  of  application  of  aj.  Tha  bara  hull  ooaffloianta 
wara  aatimatad  uaing  alandar  body  thaory  aa  davalopad  by 
Munk,*"  and  modified  by  Allan  and  Peiklna,*4  Kelly,8*  and 
Martin***  On  tha  baaia  of  thaaa  reailta, 

■  -s-° 

(«;>bh  -  2<°P  ♦  *{> 

■  2*{ 

*  -H‘ 

where  0p  la  tha  prlamatlo  ooaffioiant  and  x£  ia  tha  non- 
dlmanaionallaad  x  ooordlnata  of  tha  baaa  of  tha  hull. 

Tha  vaJue  of  waa  aatimatad  uaing  two  diffarant 
thaorlaa.  Tha  flrat  of  thaaa*  dua  to  Plax  and  Lawrenoe,** 
uaaa  low  aapaot  ratio  wing  thaory**  to  pradlot  tha  contribu¬ 
tion  of  tha  fine,  and  tha  raaulta  of  Low  and  Stone,*4 
Spreiter,**  and  Lannarta40  to  aocount  for  tha  fin  body  inter- 
faranoa  affaota.  aj  ia  datarmlnad  from 


*1  ■  o  (i  +  £)  (i  ;• 

"a  t 


Ca. 


Bi\ 


(Aj 


P 


(B-5) 


whara  C  la  tha  fin  normal  foroa  derivative  ooaffioiant 

*a 

baaad  on  wing  area,  d  la  tha  body  diamoter,  bfc  la  tha  tip-to- 
tip  fin  epan,  C_  ia  tha  body  intarfaranoa  normal  foroa  oo¬ 


affioiant,  C  la  the  fin  normal  foroa  ooaffioiant,  (A  )p  la 

**  * 

tha  projected  area  of  tha  (<?)  finu  and  Ax  la  tha  body  oroaa- 

aectional  area.  Tha  valuaa  of  both  C_  and  tha  ratio  C„  /c. 

*a  *Bi  * 


w 
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are  obtained  from  curves  presented  by  Flax  and  Lawrence.  The 

values  obtained  are  C_  *  -2.5,  C„  /c_  ■  O.38,  on  the  basis 

za  ZB1  zw 

of  which,  eq.  B-5  yields  z{  ■  -11. 7. 

The  second  method  of  predicting  z[t  used  previously 
by  Davidson  Laboratory  and  others  In  airship  and  torpedo  work, 
uses  the  Prandtl  formula  for  alllptioally  loaded  high  aspeov 
ratio  wings,18 

(A  ) 

■'  “  "  14-0AR'1  TT"  (B-6) 


where  the  augmented  aspeot  ratio,  AR',  is  defined  by 


AR« 


(B-7) 


where  0  Is  the  ohord  of  the  fins,  This  definition  of  the 
aspeot  ratio  of  the  fin  system  accounts  for  the  effects  of 
all  fin-body  interactions.  On  the  basis  of  eqe.  B-6  and  B-7, 
Z 1  -  -9.6. 

No  theoretical  or  empirical  method  for  determining 
x{  precisely  In  known.  Experimental  evidence  indicates,  how¬ 
ever,  that  x|  corresponds  to  a  longitudinal  coordinate  of  the 
fin  on  the  body.  Accordingly,  estimates  for  M^,  Z',  and 
were  made  on  the  basis  of  all  x{  values  corresponding  to  fin 
coordinates,  and  for  z[  values  derived  using  both  the  low 
aspeot  ratio  thecry  and  the  airship  theory.  The  results  are 
presented  in  graphical  form  in  Fig.  3.  The  values  ultimately 
selected  for  use  in  the  trajectory  predictions,  corresponding 
to  z{  ■  -9.6  and  x{  ■  -0.31,  also  are  listed  in  Table  I. 

The  drag  coefficient  X£,  was  estimated  with  the 
guidance  of  empirical  results  published  in  the  literature. 

The  effect  on  the  drag  coefficient  of  the  mutual  interference 
between  the  body  and  the  fins  was  assumed  negligible.  Thus, 


(B-8) 
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where  and  (x*)p  are,  respectively,  the  coefflolente 

for  the  bare  hull  and  for  the  fine,  each  evaluated  ae  though 
the  other  component  were  not  present.  A  value  of  (X')p-  -0.2' 
was  estimated  on  the  basis  of  experimental  data  for  two- 
dimensional  wedges  presented  by  Hoerner,41  The  friction  drag 
coefficient  for  the  bare  hull,  (X£  )K?  ,  was  estimated  to  be 
•0.13,  a  value  whioh  approximates  Che  Schoenherr  friction 
line  for  Reynolds  numbers  near  2x10* ,  The  bass  drag  coeffi¬ 
cient  for  the  bare  hull,  (X£  )m  ,  was  estimated  using  an 
empirioal  formula  presented  By  Hoerner, 41 

(Xob*BH  "  ‘ 


whioh  gives  (X£  )BH  •  -0.08,  It  was  assumed  that  the  bare 
b 

hull  coefficient,  (x^)^,  is  equal  to  the  sum  of  (x»  and 
(X»  )w  .  On  this  basis,  eq.  B-8  yields  X^  -  -0.45. 


V® 


The  longitudinal  virtual  mass  coefficient,  X^,  also 
was  estimated  on  the  basis  of  published  empirioal  results. 
Again,  mutual  interference  effects  were  neglected,  so  that 

*J  ■  <*J>bh  +  (XJ>F  <B-9> 

where  the  subscripts  have  the  same  connotation  as  in  the 
representation  for  X£.  Yee-Tak  Yu86  reports  experiments  in¬ 
dicating  that  the  virtual  mass  of  two  cylindrical  cups  placed 
base  to  base  is  equal  to  the  virtual  mass  of  a  disk  having  a 
radius  the  same  as  that  of  the  cylinder  plus  the  mass  of  the 
water  contained  in  the  cylinder.  On  the  basis  of  these  re¬ 
sults,  it  was  assumed  that  the  virtual  mass  for  the  Basic 
Finner  bare  hull  may  be  approximated  by  the  value  obtained 
for  the  corresponding  disk.  According  to  the  experimental 
results, 26  this  is  given  by 
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B-8 


I 


[ 

[ 

[ 

l 


(Xl)_u 

uBH 


3.41  rac 
&<■ 


whenoe,  (X(I)bh  ■  -0.11 

It  was  assunsd  that  the  longitudinal  virtual  mass 
of  one  of  the  Basic  Pinner's  wedge-shaped  fins  is  equal  to 
that  of  a  flat  plate  fin  with  thiolmess  equal  to  the  maximum 
thiokness  of  the  wedge.  This  value  was  estimated  by  means 
of  an  empirical  formula  for  the  virtual  mass  of  a  rectangular 
parallelepiped  also  developed  by  Yee-Tak  Yu,*7 


*v  -  P  ’('^a  +3.5  s.s.t//*  ( B-lo) 

where  is  the  virtual  mass,  in  grams,  p  is  the  fluid  density 
in  grams/cm9,  and  2  sx,  2  sa,  and  2tx  are  the  dimensions  of 
the  parallelepiped,  in  om,  2tx  being  the  dimension  in  the 
direction  of  the  aooeleratlon.  Since  this  formula  is  not 
dimensionally  homogeneous,  it  oannot  be  written  in  terms  of 
nondlmensional  quantities.  Consequently,  the  calculation 
must  be  made  relative  to  a  particular  sized  configuration. 

For  this  purpose,  sx  -  tx  2  cm}  sa  >0.16  om  was  chosen.* 

This  corresponds  to  a  configuration  4  cm  in  diameter  (see 
Pig.  l),  Intermediate  in  size  relative  to  the  configurations 
employed  in  the  Basic  Pinner  trajectory  experiments.  On  this 
basis,  -  1.91  grams,  whenoe  (XJ)p  -  >0.03.  Substitution 
into  eq.  B-9  gives  XI  ■  -0.14, 

To  estimate  the  normal  virtual  mass  coefficient, 

ZJ,  fin-body  interference  effects  again  were  assumed  to  be 
negligible.  Thus, 


*  The  variation  in  size  of  the  configurations  employed  in  the 
model  tests  has  very  little  effect  on  the  results  of  these 
calculations. 
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( B-ll) 


H  "  <zi>BH  +  <*i>p 

The  bare  hull  coefficient,  was  assumed  to  be  equal  to 

that  of  a  prolate  spheroid  with  the  same  length- diameter 
ratio  as  the  Basie  Pinner  hull.  According  to  the  results  of 
Lamb,18  this  is  given  by 

( Z I )  Maggi  ,-0.96  2C 

w  Hi  |  At  P 

whence,  (ZJ)^  -  -1.6 

It  was  assumed  that  the  normal  virtual  mass  of  one  of  the 
Basic  Firmer 's  wedge-shaped  fins  is  equal  to  that  of  a  flat 
plate  with  thickness  equal  to  the  average  thlokness  of  the 
wedge.  This  value  was  estimated  by  applying  eq.  B-10,  in  this 
instance  taking  s1  >  sa  -  2  cm,  tx  ■  0.08  cm.*  Thus  my  -  39.6 
grams,  whence  {Z£)p  ■  -0,3.  Substitution  into  eq.  B-ll  gives 

ZJ  -  *1‘9. 

The  virtual  moment  of  Inertia  coefficient,  M^,  was 
estimated  by.  meane  of  the  equation 

+  X*‘Ui>F  (B-l2) 

Thie  equation  was  derived  by  analogy  with  eq.  B-l,  Inter¬ 
ference  effects,  however,  were  neglected.  Hence  is  the 
dimensionless  longitudinal  coordinate  of  the  point  of  applica¬ 
tion  of  (Z£)p ,  taken  to  be  the  longitudinal  midpoint  of  the 
fin,  i.e.,  x£  -  -0.35.  The  bare  hull  coefficient, 
was  assumed  to  be  equal  to  that  of  a  prolate  spheroid  with 
the  same  length-diameter  ratio  as  the  Basic  Firmer  hull. 


*  See  note  on  preceding  page. 
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From  Lamb i 


<»{> 


IB 


BH 


0.07  I, 

fv ' 


where  X  la  tha  moment  of  inertia  of  the  displaced  fluid. 
yo 


Thus 


Mi  - 


•0.08.  Substitution  into  eq.  B-12  gives 


The  eecond  order  derivative  X*  was  estimated  on 

wq 

the  basis  of  the  results  of  potential  theory*  whioh  indicate1" 


X»  -  ZJ 
wq  w 


whence* 

x;q  - 

The  values  of  each  of  the  hydrodynamlo  coefficients 
derived  hereiu  are  reproduced  in  Table  V. 
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APPENDIX  C 

SOLUTION  OP  VERTICAL  PLANAR  MOTION  EQUATIONS 

Numerical  (Exact)  Solution 

Let  tQ  denote  the  initial  value  of  the  independent 
variable  t  ( time) .  The  initial  valuee  u( tQ) ,  w( tQ) ,  q( tQ) , 
and  ®(tQ)  are  given.  Equations  2  through  4  are  solved  for 
the  three  unknown  quantities,  u(tQ),  w(tQ),  q(tQ),  to  obtain 

4<t )  -  ([“<«■>>] *♦[»<»»>]•>•  i -  gSffjJ «<»,>«(»,>  ♦  ■»  [«<»„>]*  t  n- 1  \ 

*  Ml 


[«*0>]  ♦  ♦  [-  ♦  K«j]  *  fv*"0 

•’ - — - — ' 

.  '{"* * '  K**yiu(y  «<o  •  % «“  t*(t0» } 

_  [*  *  ^a*][I»»  ’  "V  (0-2) 


4(U 


^^(Q^Q  -[wq-  Kt,Mo.  “(Q^O  -  ^  0011  ,*(Q.  +  ”o*<0 


yy 


(C-3) 


where  the  value  of  w( tQ)  ie  first  obtained  from  eq.  C-2  and 
then  eubstituted  into  eq.  C-3  to  solve  for  q(tQ).  The  values 
of  the  linear  and  angular  velocities  at  time  tQ  +  At  are,  to 
the  first  order  of  approximation, 

u(tQ  +  At)  -  u(tQ)  +  u(tQ)At  (c-4) 

w(tQ  +  At)  -  w(tQ)  +  w(t0)Afc  (C-5) 

q(tQ  +  At)  ■  q(tQ)  +  q(tQ)At  (C-6) 
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Now,  Binoe  0  -  0,  q  *  9  ( bob  Appendix  A) .  Thus 

&(  fcQ  +  At)  -  e(t0)  -f  &(tQ)At  -  e(t0)  +  q(t0)At  (c-7) 

The  veloeltleB  of  the  body  relative  to  the  fixed 
(xQ,  yQ,  sQ)  coordinate  systems  are  obtained  from  the  results 
of  eqs.  C-4,  C-5,  C-7  by  applying  the  following  transforma¬ 
tion  equations  (see  Table  of  Direction  Cosines,  Appendix  A)t 


*0(t)  -  u(t)  cos  8(t)  +  w(t)  sin  0(t) 

(C-8) 

iQ(t)  -  -u(t)  sin  fl(t)  +  w(t)  cos  0(  t) 

(C-9) 

The  position  is  then  obtained  from 
x0(tp  +  At)  -  x0(tQ)  +  x0(t0)At 

( c-10) 

z0(t0  +  At)  -  a0(t0)  +  s0(to)At 

(c-ll) 

where  the  initial  values  xQ(to)  and  zQ( tQ)  are  given  quanti¬ 
ties. 


The  results  of  eqs.  C-4  through  C-ll  are  now  re¬ 
garded  as  the  initial  data,  and  the  same  procedure  is  re¬ 
peated.  The  process  is  continued  until  t  attains  a  value  cor¬ 
responding  to  the  value  of  time  specified  as  the  endpoint  of 
the  run. 

If  the  functions  involved  are  well  behaved,  the 
sequences  of  functions  obtained  by  this  procedure  will  con¬ 
verge  as  the  increment  At  1b  taken  smaller  and  smaller.  In 
the  present  case,  the  reaultB  obtained  on  the  baals  of 
At  »  0.005  sec  are  indistinguishable  from  those  computed 
using  At  =  0,0025  sec.  Accordingly,  the  former  value  is  used 
in  the  predictions. 

The  solution  was  programned  and  coded  using  the  IBM 
1620  Fortran  programming  language.37  The  calculations  were 
performed  using  the  IBM  1620  Computer  Stevens  Institute  of 
Technology. 
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Analytic  (Approximate)  Solution 

Upon  substitution  of  eqs.  9  and  10,  eqs.  2  through 
4  may  be  written  In  the  form 


"  -  fvxj 

w  -  A  +  (Bw  +  Cq)u 
q  -  A*  +  (B*w  +  C*q)u 
where, 

( W-B)  ( 3!yy  -  §AxlaMl)  -  Wmx0a 
A - B - - - 

r  &?i  -  fv3||j  +  fvs  mo 


(C-12) 

(C-13) 

(C-14) 

(C-15), 

(C-16) 


("*  f*xl2i>(ljry  •  -  "V™0  . 

C  ■  . . .  . . . .  . . b . . . . . . . . . . .  IC-17/ 

(W-B) mxQ  -  WxQ(m  -  |AxtZl) 

A*-  — ■ - jj - H-  (C-18) 

f‘xzi"”‘<l  +  fV“i< "  -  £*,“<>  ,  , 

B*- - £ - —  (C-19) 

(nrt*  |Ax-tZ‘)mx0  +  (fAxt8M*  -  mxQ)(m  -  |AxtZ!) 

c*« - - - p— 2 ~  ( c-2o) 


D  -  (m  -  fVZi)(lyy  '  “  "^Q  (C-2l) 

The  longitudinal  force  equation  (eq.  C-12)  is  identical  to  the 
equation  for  the  straight  horizontal  trajectory  (eq#  6), 

Hence  the  solution  is  given  by  eq.  11.  This  is  substituted 
into  eqs.  C-13  and  C-14.  Introduction  of  a  new  independent 
variable. 
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r  ■  bt  +  1 


(C-22) 


in  the  resulting  equations  gives 

b  -  A  +  £  (Aw  +  Mq)  (C-23) 

b  -  A#+  i  (x\  +  n*q)  ( C-2 4) 

where 

X  -  Bu(0),  n  -  Cu(o) ,  \*  m  B*u(  0) ,  fi*  -  C*u(  0)  (C-25) 

Equation  C-24  is  solved  for  w  and  differentiated  with  respeot 
to  r.  The  results  are  substituted  into  eq.  C-23  to  obtain 

b*T*  0  +( b-M*-x) br  gSL  +( q-( A*A+bA*-AA*) r  ( C-26) 

The  homogeneous  equation,  a  differential  equation 
of  the  Euler- Cauchy  type,  has  the  solution, 

°x  o8 

qH  -  KA  t  +  K,j  t  (C-27) 

where 


„l(,  .  (X  *  M«)  tylu  -  4(„.>  .  ux.)  (0-28) 

The  particular  solution  is  ^iven  by 

qp  "  *8  r  (c-29) 

where 

„  _  X*A  +  (b  -  X)A*  , _ . 

Ks  b=---rr+  V)  +  Vim  -  mx*)  (c'-o! 

Hence  -he  complete  solution  is 
ox  a3 

q  =  Kx  r  +  K2t  +  KgT  (c-3l) 

as  can  be  verified  directly  by  substitution  into  eq,  C-26. 
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Equation  C-31  la  substituted  into  aq.  C-24  to  ob¬ 
tain  tha  aolution  for  tha  normal  valooity, 

•  -  RH  </ *♦  m  & 

(C-32) 

The  ramaining  unknown  constants,  Kx  and  Kg,  detar- 
mined  from  the  initial  oondltiona,  are  given  by 


K,  - 


i  r  i  * 

l-°i  b°i-u#  ,  .  .* 

r-  Ka+  _ _ SL  q(o)  -  -y  -£1L-  w(0)-r— A- 

°i"°a  fe(ox-ot)  b(o4-  era  bToj-c 


Kx  -  q(  0)  -  Ka  -  Kg 


(C-33) 

(C-34) 


Jd  JkQ 

Now,  ainoe  0-0,  q  «  ^  -  b  g~.  Therefore,  integrating 
eq.  C-31, 


9  m 


Eto7tir]T°1+1  +  [ftIJtd]  t°"+1  ♦  h  +  +  *.  <°-35> 


whei'e  the  constant  of  integration,  determined  from  the  initial 
conditiona,  la  given  by 


K4  -  0(0) 


(C-36) 


Solutiona  for  the  longitudinal  and  vertical  dis¬ 
placements  may  be  obtained  similarly  by  integrating  eqa.  11 
and  C-32.  This  is  not  done  in  the  present  work,  alnce  it  is 
fait  that  the  solutiona  already  obtained  are  sufficient  for 
the  comparisons  for  which  they  are  Intended.  The  important 
results  are  repeated  below  in  terms  of  the  original  independ 
ent  variable,  t. 


u 


u(0) 

birr 


(n) 
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*  ■  (bOi-i^KiOxs+l)01'1  +  (bo,-  fi*)  JC,(bt+l)9,**+  (b-(i«)K,-A*  (bt+1) 

(C-37) 

6  -[l^]<b*«>Bl%tikl)]  (»t«)"*+1+  &MH)'  *  K. 

(C-38) 

The  various  constants  appearing  in  the  equations  are  evaluated 
by  the  eucoessive  application  of  the  equations  by  which  they 
are  defined.  The  velooJ  cies  of  the  body  relative  to  the 
fixed  (xQ,  yQ,  zQ)  coordinate  system  are  obtained  from  the 
above  results  by  applying  eqs.  C-8  and  C-9. 
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TAILS  Ill 

COMPARISON  OF  VALUES  OF  z',  PREDICTED  USING  TWO  DIFFERENT  THEORIES 
WITH  EXPERIMENTALLY  DERIVEO  VALUES 

CASE  NO. 

VEHICLE 

REFERENCE 

VALUES  OBTAINEO  FOR  z', 

EXPERIMENTAL 

FLAX  AND  LAWRENCE 
THEORY 

■ 

1 

Submarine 

38 

-2.78 

-3.67 

-3.15 

2 

Misti  la 

29 

-2.29 

-3.24 

-2.51 

3 

Misti  la 

39 

-2.51 

-2.01 

-2.44 

4 

Misti  la 

22 

-14.5 

,  —  .  - . 

-13.5 

-11.8 
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TABLE  V 

PARTICULARS  OF  THE  BASIC  FINNBR  MISSILE 


Dlanatar  o f  bar*  hull,  d  1.000 

Langth  of  bar*  hull,  I  10.00  d 

2 

Maxlaua  croaa-aaetlonal  area  of  bar*  hull,  Ax  0.785  d 

Nattad  araa  of  bar*  hull,  S^  27.0  d2 

Baa*  araa  of  flna  (4),  (Afe)p  0.320  d2 

Nattad  araa  of  flna  (4),  Sp  8.10  d2 

Projaotad  araa  of  flna  (2),  (Ap)p  2.00  d2 

Total  voluna,  Vol.  6.33  d3 

Moaant  of  lnartla  of  fluid  dlaplaead  by  hull,  I  30.8  pd' 


Hydrodynaalc  Coafflctnntai 


X'  . 0.43 

0 


xi 


-0.14 


W-  0 


V  -V*  *1<9 

VXrr-  0 

■  YJ . 11.0 


•  Y!....  -1.9 


Z'  -  -Y'...  -3.78 
d  * 

ZJ.-Yf...  0 

-  -NJ...  -2,12 


M^  -  -N|...  0 


M4 


NJ . 1.24 

NJ....  -0.12 


*  So*  aaction*  baginning  on  pp,  10  and  19, 


TABLE  VII 

COMPARISON  OF  DL  ANG  NFS  RESULTS 


H 

CIT 

RUN 

NO. 

MAXIMUM  DIFFERENCE  BETWEEN  OBSERVED  ANO  COMPUTED  VALUE 

LONOITUOE,  Ktf 
MOOEL  LENGTH 

OEPTH,  «# 

MODEL  LENGTHS 

INCLINATION,  B 

OEGREES 

OL 

OL 

NPO 

OL 

NPG 

a 

F- 96 

.09 

.11 

.07 

.ia 

1.8 

m 

8 

F^a 

.11 

.19 

.01 

.06 

0.4 

i.i 

10 

F- 64 

.07 

.00 

.01 

.08 

0.4 

0.8 
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FIGURE  L  BASIC  F1NMER  MISSILE 
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REYNOLDS  NUMBER  X  10 


FIGURE  2.  EXPERIMENTALLY  DETERMINED  VALUES  OF  DRAG 
COEFFICIENT  VERSUS  REYNOLDS  NUMBER 
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tflSl,  DECREES 


O  t\ ,  MSEC  ON  AIRSHIP  THEORY 
4  l! ,  BASED  ON  tOW  ASPECT  RATIO  WINS  THEORY 
—  NPO  CURVE-PIT" 


} 


SEE  TABLES  POR  VALUES  OP  HYOROOVNANNC 
COEFFICIENTS 


CIT  RUN  P-SS 
(OL  CALCULATION  I) 


AVERAGE  POR  THREE  RUNS 
CIT  RUNS  F*SG,8B,G4 
(DL  CALCULATIONS  *,S,  10) 


S 

o 

X 

3 

d 


0.3  - 

0.2  - 
0.1  " 


J> 

4 


2. OS 


AVERAGE  MEASURED  COEFFICIENTS 
|  SEE  TABLE  S 

4 

- 1.  — +—  —  —  JX- 
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AVERAGE  MEASURED  COEFFICIENTS 
SEE  TABLE  H 
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FIGURE  4.:  MAGNITUDES  OF  MAXIMUM  DIFFERENCES  BETWEEN  OBSERVED  AND 
COMPUTED  VALUES  OF  TRAJECTORY  PARAMETERS  WITH  VARIOUS 
SETS  OF  HYDRODYNAMIC  INPUT  DATA 
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FIGURE 


5.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  TRAJECTORIES  FOR 
QL  CALCULATION  NO.  I 


p-flq8 


DEPTH- INCHES,  (-HNCUNATKM-OEGftEES 


VELOCITY  -  EPS 


FIGURE  6.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  TRAJECTORIES  FOR 
D.L,  CALCULATION  NO  2 
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DEPTH- INCHES.  (-)  MCLMATKM.  DEGREES 
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FIGURE  7,  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  TRAJECTORIES  FOR 
D.L.  CALCULATION  NO,  3 
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OEPTH-WCHES.  (-»  MCUNATKM.  DEGREES 


FIGURE  9.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  TRAJECTORIES  FOR 
DL.  CALCULATION  NO.  5 
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PTH- INCHES,  (-)  INCLINATION.  bEOKEES 


FIGURE  10.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  TRAJECTORIES  FOR 
D.L  CALCULATION  NO,  6 
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DEPTH- INCHES.  (-)  INCLINATION.  DEGREES 
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FIGURE  I!.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  TRAJECTORIES  FOR 
D.L.  CALCULATION  NO,  7 


(-)  depth,  inches,  inclination,  oeonei 


FIGURE  12.  COMPARISON  OF  EXPERIMENTAL  ANO  PREDICTED  TRAJECTORIES  FOR 
D.L.  CALCULATION  NO.  8 
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(-)  DEPTH,  INCHES.  INCLINATION,  DEGREES 
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FIGURE  13.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  TRAJECTORIES  FOR  D.L.  CALCULATION 
NO.  9 
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FIGURE  15.  COMPARISON  OF  EXACT  (SOLID  LINES)  AND  APPROXIMATE 

(DASHED  LINES)  SOLUTIONS  FOR  FOUR  TYPICAL  HIGH  SPEED 
VERTICAL  TRAJECTORY  PREDICTIONS. 
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FIGURE  16  COMPARISON  OF  EXPERIMENTAL  ANO  PREDICTED  TRAJECTORIES  FOR  D.L  CALCULATION 
NO.  19 
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FIGURE  17.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  TRAJECTORIES  FOR 
D.L.  CALCULATION  NO-  12 
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FIGURE  18.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  TRAJECTORIES  FOR 
DL  CALCULATION  NO.  16 
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FIGURE  19.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  TRAJECTORIES  FOR 
n«_.  CALCULATION  NO.  18 
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FIGURE  20.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  TRAJECTOf  cS  FOR 
DL  CALCULATION  NO  20 
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FIGURE  21.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  TRAJECTORIES  FOR 
D.L.  CALCULATIONS  NOS.  21,  22,  23 
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FIGURE  2  2,  COMPARISON  OF  EXPERIMENTAL  ANO  PREDICTED  TRAJECTORIES  FOR  DL 
CALCULATIONS  NOS  24,  25,  26 
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